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ABSTRACT 


Mean velocity profile data are reported for blown, 
unblown, and sucked accelerated boundary layers. The pressure 
gradients investigated are those corresponding to constant 

V dU oo 

values of the pressure gradient parameter K = — tt -t— . The 

Uco 

three values of K considered are 0.57 x 10 - ^ , 0.77 x 10“^ , 
and 1.45 x 10”^ . For each pressure gradient, the surface 
boundary conditions c /er a range of constant blowing and 
sucking fractions from F = -0.004 to 0.006 . 

The friction factor was found to be decreased by the 
imposed acceleration for suction, whereas for blowing an 
increase was indicated. This is based upon a comparison 
with zero pressure gradient data for equivalent local momentum 
thickness Reynolds number and blowing parameter B . 

Velocity profiles corresponding to these accelerated 
flows were found to deviate from those characteristic of 
zero pressure gradient flows. These deviations in the fully 
turbulent and wake regions of the boundary layer apply over 
the range of blowing and sucking fractions investigated. 

For boundary layers where local momentum thickness 
Reynolds number is invariant in the flow direction, velocity 
distributions near the wall were found to be similar in 
U + vs y + coordinates. In the exceptional cases where 
similarity was not observed, substantial structural changes 
in this region are suggested. 

The corresponding velocity distributions far from the 
wall were found to attain similarity in U/U^ vs y/6 coor- 
dinates and velocity-defect coordinates. These outer region 
coordinate systems are equivalent in these cases, as a result 
of the relative constancy of friction factor in the flow 
direction. 
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Unique values of local momentum thickness Reynolds 
number, Regg » and shape factor, H , were found to be 
associated with these boundary layers for given value of the 
parameter K and blowing or sucking fraction P . 

Semi-empirical representations of these data in the form 
of simple Prandtl mixing length distributions are presented. 
The inner regions of the layer are represented by means of a 
two-layer model and a continuous model. The outer regions 
are represented by an appropriate truncation of these dis- 
tributions. 

Utilizing the resulting eddy-diffusivity relations in 
a finite-difference procedure, predictions of the present 
data were obtained. Agreement with the data is found to 
be acceptable for most engineering purposes. 
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CHAPTER I 


INTRODUCTION 


The turbulent boundary layer with non-zero normal 
velocity V w at the surface is of practical interest since 
injection and suction of fluid at a surface are utilized for 
thermal protection and boundary layer control, respectively. 

In most of these applications, the mainstream fluid is locally 
accelerated or decelerated; hence, the combined effects of 
transpiration at the surface and acceleration of the main 
flow must be considered in system design. 

Since t' ere is no adequate tneory for shear flow turbu- 
lence, existing turbulent boundary layer "theory" relies 
heavily on experimental results. Well documented velocity 
measurements are necessary to extend the theory so that these 
combined effects are properly represented. 

The present investigation is concerned with the fluid 
dynamics of accelerated turbulent boundary layers with in- 
jection or suction of fluid at the surface. Experiments are 
restricted to boundary layers characterized by constant 


values of the acceleration parameter KlK - 


V dU oo 


blowing fraction FjF = 


te) • 


U„ dx 


) 


and 


The specific flows con- 


sidered are further characterized, as: two-dimensional, 

constant property flow over a uniformly permeable surface, 
where the pore opening and pore spacings are small enough to 
insure an aerodynamic ally smooth surface with uniform in- 
jection or suction. 


A. Review of Previous Experimental Works 

In recent years a number of experimental hydrodynamic 
investigations have been concerned with blown (V > 0) and 
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sucked (V < 0) turbulent boundary layers [1-17]. Only two 
of these investigations are known to consider the effects of 
acceleration of the main flow. Experimental hydrodynamic 
investigations of accelerated boundary layers on impermeable 
surfaces have also been reported [18-24]. Results of these 
experiments relevant to the present study are briefly dis- 
cussed. 


A.l. Cases of constant free-stream velocity with 
transpiration 

Simpson [17] presents a detailed review of the 
experimental investi 3 ations, up to 1967 , in the case of zero 
pressure gradient. The reader is referred to this review 
which will not be repeated here. 

The experimental data of Simpson, obtained on the same 
apparatus as used for the present investigation, are believed 
to represent a reliable and consistent set of hydrodynamic 
data covering the widest range of both blowing and sucking 
conditions. These data are also documented in a manner suf- 
ficient for purposes of comparison and correlation with the 
present data. As a result, these data and associated cor- 
relations proposed by Simpson are utilized in the present 
study. 

A. 2. Accelerating flows without transpiration 

The relevant experimental investigations of ac- 
celerated flows on impermeable surfaces may be categorized 
as: l) mild pressure gradient flows of the turbulent 

"equilibrium" type, or 2) strong pressure gradient flows 
where the onset of re laminar izat ion is of major concern. 

The experiments of Herring and Norbury [18] represent 
one of the major contributions in mild pressure gradient 
turbulent boundary layer investigations. Velocity-defect 
distributions for accelerated boundary layers were established 
corresponding to values of the dimensionless pressure gradient 
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parameter (3(£ - b,/T dP/dx) of -0,35 and -0.53 . These 
velocity-defect profiles are shown to be similar in tne outer 
regions of the layer, as previously proposed by Clauser [ 25 ]. 
Their data agreed well with the theoretical predictions of 
Mellor and Gibson [26], which are based upon "equilibrium" 
theory. It is shown in section B of this chapter that the 
unblown asymptotic boundary layer flows considered m tne 
present study are also constant flows where £ < - 0,5 . 
These findings suggest that "equilibrium" theory might pos- 
sibly extend to the strongly accelerated flows presently 
considered . 

The experiments of Launder et al. [19-21] nave been con- 
centrated in tne area of strongly accelerated boundary layer 
flows. The major concern in these investigations has been 
to determine the onset of relaminarization. The majority of 
the investigations of Launder et al. were concerned with 
constant K accelerations, shown to be an important parameter 
by several studies. The thermal measurements made by Moretti 
and Kays [ 27 ] for strong, constant K , accelerated flows 
indicate very substantial decreases in Stanton number wnich 
could be roughly correlated with the parameter K , In con- 
current studies made by Kline, Reynolds, Schraub, and 
Runstadler [22] it was demonstrated, that the burst frequency 
of turbulent disturbances originating at the wall also cor- 
related with the parameter K . 

Launder and Stinchcombe [20] studied constant K ac- 
celerations where tne local momentum thickness Reynolds 
number, Reg 0 , approached an asymptc tic limit. It was 
observed that the velocity profiles were similar in such a 
boundary layer for values of K = 0.7 x 10 - ^ , 1.25 x 10 - ^ 
and 3 x 10~^ . Velocity profile similarity is predicted 
by the exact laminar solution* for constant K accelerations. 


The reader is referred to reference [28] for a statement of 
indicated laminar solution. 
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suggesting similarity through values of K corresponding to 
complete relaminarization. Profiles for these strongly ac- 
celerated flows are shown to differ from those corresponding 
to turbulent boundary layer flow over an impermeable flat 
plate. As K was increased, a continuous shift from a 
typical turbulent profile was exhibited. This shift is char- 
acterized by an "overshoot" of the velocity profile from that 
corresponding to the "law of the wall" and a simultaneous de- 
crease of the wake region. Later experiments of Launder and 
Jones [21] qo not corroborate the quantitative results of 
Launder and Stinchcombe even though performed on the same 
apparatus, but the same qualitative conclusions were found 
to apply. This profile "overshoot" is consistent with the 
findings, presented in reference [22], that a decrease in 
the burst rate of turbulent disturbances originating at the 
wall is associated with an increase in K , resulting in a 
reduction of turbulent fluctuations at the wall. 

The experiments of Badri Naroyanan and Ramjie [23], 
concerned with both constant and variable K flows, demonstrate 
the same profile "overshoot" behavior. 

The experiments of Patel and Head [24] were concerned 
with boundary layer flows for which K was strongly varying. 

The characteristics observed in the initial adjustment of 
the layer to the imposed pressure gradient are noteworthy. 

For accelerations where the parameter K is increasing in 
the flow direction, the profiles are shown to initially 
"undershoot" the "law of the wall" and then approach the 
previously observed "overshot" appearance in a continuous 
development. Although corroborative experimental data is 
not known to exist, this behavior is thought to be indicative 
of the inability of similarity data to adequately represent 
boundary layer flows in which K is strongly varying. Non- 
similar flows of this type are of importance, but are not 
treated in the present study. 
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A . 3 . Accelerating f 7 ows with transpiration 


The characteristics found in accelerated boundary- 
layers on impermeable surfaces will undoubtedly be present 
when blowing or suction at the surface is introduced. In the 
work of Simpson [ 17 ] > the effects of transpiration in zero 
pressure gradient flows are shown to be quite substantial in 
terms of order of magnitude reduction in friction factor 
with blowing and corresponding increases with suction. Mean 
profile appearances are also shown to significantly deviate 
from those of flow over an impermeable flat plate; upward 
shifts of the velocity profile from the "law of the wall' 1 with 
blowing and downward shifts with suction are indicated. 

Although these characteristics a.-sociated with transpiration 
are also expected in the combined case, any successful attempt 
to represent the boundary layer characteristics would be only 
fortuitous if based on data where each effect has been isolated. 

Only two experimental works are known to exist v/here tne 
combined case was investigated. Each considered only blown 
layers and an insufficient amount of experimental data to 
adequately represent the boundary layer characteristics. 

Romanenko and Kharchenko [7] present friction lactor and. 
Stanton number data but no profiles of velocity or temperature. 
Insufficient documentation of the conditions under which these 
data were obtained in the combined cases, including a speci- 
fication of the free-stream velocity distribution, prevents 
its utilization in the present study. 

In the experiments of McQuaid [10], friction f .ctors 
were not measured in the two combined biowing and acceleration 
runs made. This presents a restriction on the formulation of 
data correlations. Utilizing friction factors corresponding 
to Stevenson's inner law [29], McQuaid was able to predict 
momentum thickness distributions which agreed well with the 
experimentally determined distributions. Although this agree- 
ment is found to exist, this is considered not to be a 
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sensitive test of friction factor when blown boundary layers 
are considered. 

It is concluded that a need still exists for a complete 
set of documented data for the combined case, where friction 
factors are accurately known. 


B . Description of Asymptotic Boundary Layer Flow 

The asymptotic boundary layer flows considered in the 
present study are charactei ized by acceleration at a constant 
value of K and transpiration at a constant value of F . 

To illustrate the necessity of these boundary conditions, the 
two-dimensional momentum integral equation can be presented 
in the form 


dRe*, 

dR 


= 2=- - Reg (H+l)K + F 


(I-D 


where dR 


x 


U. 

- 


For constant values of K and F , the boundary layer develops 
such that the terms on the right-hand side of equation (I-i) 

dRs§2 


balance, forcing the derivative 


dR. 


to zero. Such a 


boundary layer will be classified as asymptotic in the regime 
where Res^ is maincained constant. Note that this condition 
requires a constant friction factor C for similar profile 
development . 

These asymptotic flows are also constant f3 flows in 
the event similar profiles exist. For zero blowing, sub- 
stitution of the asymptotic form of equation (I-l) into the 
definition of £ vields 


P = - 


H 

H+l * 


( 1 - 2 ) 
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The data of Launder and Stincncombe [21] show the existence 
of similar profiles for asymptotic accelerated boundary layers 
where equation (1-2) indicates that tnese are constant £3 
flows where £3 < -0,5 - 

The present study is restricted to these asymptotic 
boundary layers for purposes of convenience. For asymptotic 
or near-asymptotic flows, equation (I-l) yields one method 
of estimating friction factor, where the derivative, 
dRe&p 

-T 5 , represents a correction to the asymptotic form of 

ok x 

equation (I-l) - This is a desirable characteristic since tne 
direct measurement of friction factor is not available on the 
present apparatus. When the friction factor remains constant 
these flows are also characterized by constant values of the 
blowing parameters B and V 4 " , as well as p + , These 
characteristics are desirable in the formulation of data 
correlations 

C . Objectives of the Present Research 

The motivations for the present research can be summarized 
as follows. Turbulent boundary layer theory is heavily depen- 
dent upon the results of experimental studies where specific 
effects have been isolated. Several experimental investiga- 
tions have been reported where the effects of transpiration 
and the effects of strong acceleration of the main flow are 
considered independently. The few experimental investiga- 
tions where the combined effects have been studied are in- 
sufficient to represent the boundary layer characteristics 
resulting from the coupling of these effects. 

The overall objective of the present work is to investi- 
gate and record the fl\iid dynamic behavior of the turbulent 
boundary layer where the combined effects of transpiration 
and acceleration are present. The range of blowing, suction, 
and acceleration considered is that covering most practical 
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applications where turbulent boundary layer theory is 
appropriate. The following subdivisions of this objective 
are : 

1. To tabulate and document mean velocity profile data 
taken on the Stanford Heat and Mass Transfer 
Apparatus . 

2. To obtain skin friction results from the mean velocity 
profiles and estimate the reliability of these data. 

3. To examine the mean velocity profiles for develop- 
ment and similarity characteristics in asymptotic 
boundary layer flows. 

4. To obtain semi-empirical representations of these 
results in the form of simple "equilibrium" eddy- 
diffusivity models. 

5. To incorporate the resulting eddy-diffusivity rela- 
tions into a finite-difference procedure that will 
analytically reproduce the experimental data and 
provide a basis for prediction of situations not 
covered by the experiments. 

In the chapters to follow, the fluid dynamic character- 
istics of the present apparatus are briefly described as well 
as the velocity profile instrumentation. Then the experimental 
velocity profiles are presented, and data necessary to de- 
scribe the conditions under which they were obtained are 
documented. Finally, the indicated eddy-diffusivity rela- 
tions and the resulting predictions of the present experi- 
mental data are presented. 
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CHAPTER II 


EXPERIMENTAL APPARATUS 

A. Brief Description of the Apparatus* 

The apparatus used in the experiments consists of a 24- 
segment porous plate, 8 feet long and 18 inch wide. The 
plate forms the lower surface of a test duct of rectangular 
cross-section, 20 inches wide and 6 inches high at the inlet 
end of the duct. The upper surface is adjustable to achieve 
the desired velocity distribution along the duct. The 1/4- 
inch thick plates are smooth to the touch and are uniform in 
porosity within +6 percent in the six inch span centered on 
the test duct centerline, where velocity profiles are taken. 
Separate mainstream and transpiration blowers provide the 
system with air, while heat exchangers are used to control 
air temperature. 

Conventional temperature and flow rate instruments 
were used to control the operation of the apparatus. Mean 
velocity profiles were taken with stagnation pressure probes 
and manual traversing equipment. The probes used for boundary 
layer surveys have a flattened mouth 0.012 inch by 0.035 inch. 
The dynamic pressures were measured with calibrated inclined 
manometers. The probes were attached to traversing instru- 
ments fastened to a rigid support frame. These spring-loaded 
micrometer-driven instruments provide for the change and 
measurement of the probe distance from the test wall. 

B. General Physical Arrangement 

The Heat and Mass Transfer Apparatus is a two-story 
facility with the operating controls and heavy hardware 


A complete description of this apparatus is contained in 
references 37, 30, and 31. 
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on the ground level, and the test section on the second floor 
of a 15-foot tower. The first floor operating area is ap- 
proximately 10 by 12 feet, as is the deck comprising the second 
floor. Photographs of the operating console area, and of the 
test section area, are shown in Figures 1 and 2, respectively. 

The objective of the apparatus is to provide and control 
three general systems: main air stream over the test plate, 

transpiration air through the test plate, and electric heater 
power to the test plates. This is accomplished by the system 
shown schematically in Figure 3> the Test Apparatus Schematic. 

Those characteristics of the air systems relevant to the 
reliability of the present experimental data are briefly 
discussed . 


B.l. Main air system 

The flow path of the main air stream is as follows: 
(l) inlet air filter, (2) mainstream flow control valve, 

(3) main blower, (4) mainstream heat exchanger, (5) screens, 
(6) plenum chamber and primary nozzle, (7) test section. 

This system provides a velocity of approximately 44 ft/sec 
at the primary nozzle exit at full flow. 

The heat exchanger located at the upstream end of the 
nozzle acts as a flow straightener, as well as a means of 
control on the air temperature. Measurements taken by 
Moffat [30] in the entrance end of the test section show that 
the temperature is uniform within 0.25 °F, under usual oper- 
ating conditions, in the core of the mainstream. 

The plenum chamber and primary nozzle provide a two- 
dimensional contraction from the heat exchanger to the test 
section. The area ratio is 4:1, and the nozzle is designed 
for uniform acr~leration of the flow to yield a uniform 
velocity profile at the entrance to the test section. 

Boundary layer suction is not used on this apparatus. 

Initially, a 70-mesh brass screen covered the flow area 
downstream of the heat exchanger. It was later found 
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necessary to replace this screen with a double screen system, 
to attain a more uniform boundary layer development within 
the nozzle for the strongest pressure gradient flows investi- 
gated. The double screen system consists of a single 100- 
mesh and a single 50-mesh screen separated 3/4 inch from one 
another. 

A trip, of 3/8-inch wide 50(1) coarse grid carborundum 
garnet paper, is located at the exit edge of the transition 
section joining the nozzle to the test section. The trailing 
edge of the trip is located 1/8 inch upstream of the porous 
test section. When the double screen system was installed, 
an additional trip was fastened to the lower wall of the 
nozzle 17 inches upstream of the porous test section. The 
new trip consists of a l/l6-inch high by 1/4-inch wide 
bakelite trip and it was positioned to aid in the development 
of a uniform turbulent boundary layer at the trailing edge 
of the downstream trip, with a minimal momentum thickness. 
Each trip is extended the full width of the duct. 

The test section consists of three plexiglass pieces: 
two side walls, fastened to the test plate structure, and an 
adjustable upper surface. Different top covers are available 
to accommodate different free-stream velocity and blowing (or 
sucking) distributions. 

B.2, Transpiration system 

The flow path of the transpiration system is as 
follows: (l) inlet air filter, (2) transpiration blower, 

(3) transpiration flow heat exchanger, (4) transpiration 
flow header, (5) flow meter bank, (6) delivery tubes to the 
individual plate sections. 

The transpiration heat exchanger allows control of the 
transpiration air temperature, Operating the entire tran- 
spiration header system at a uniform temperature equal to 
the ambient temperature prevents heat transfer between the 
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transpiration system and the ambient environment. A single 
temperature measurement in the entrance of the header is then 
sufficient to determine the temperature of the flow through 
all twenty-four flowmeter systems fed from the header. Traverses 
across the delivery duct, at the entrance to the header, have 
shown that the temperature there is uniform within 0.3 °Fj 
consequently, a single point temperature measurement is 
representative of the flow condition. 

The rig can be converted from blowing to suction by a 
reversal of the twenty-four connecting tubes and blower con- 
nection. There is no temperature control in the suction mode, 
but none is needed due to the high heat exchanger effectiveness 
of the riser tubes. This results in all suction air being 
brought essentially to room temperature before it reaches the 
flowmeters . 

All tube and hose connections were carefully inspected 
for leakage to insure measurement of flow through each plate 
was that through the associated flowmeter. 

C . The Porous Plates and Test Plate Assembly 

The porous plates used in this apparatus are made of 
sintered bronze material. They have the following character- 
istics : 

1. Overall dimensions: 18.0 x 3.975 x 0.25 inches 

2. Porosity: approximately 40 percent 

3. Porosity and flow rate uniformity: within +6 percent 

in the center 6-inch span 

4. Surface roughness: maximum of 200 microinches (RMS) 

measured with a stylus of radius 0.0005 inches 

5. Particle shape: spherical (estimated 99 percent of 

particles ) 
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6. Particle sizes; maximum diameter 0„007 inches, 
minimum 0,0023, average diameter 0,005 inches. 

No requirement exists that all plates have the same 
porosity, simply that each plate be uniform across the center 
6-inch span. In actual use, the transpiration flow rate is 
individually measured for each plate in the apparatus, and 
ample pressure is provided by the transpiration blower to 
take care of porosity differences from one plate to another. 

The test plate assembly forms the bottom surface of the 
test section. It is 20 inches wide and 96 inches long in the 
flow direction. It is made of four subassemblies bolted to 
a common support structure. Each subassembly consists of an 
aluminum casting accommodating 6 individual porous plates. 

A cross-section view of a typical compartment in the 
test plate assembly is given in Figure The use of pre- 
plates and honeycomb on the underside of each plate was 
found necessary to insure a uniform temperature in the air 
flowing through each plate. Five thermocouples are embedded 
in each plate to measure the surface temperature. 

Porosity mappings of each plate were obtained by 
Moffat [30] o These mappings are used to calculate that 
fraction of the flow which passed through the center span of 
each plate. Small corrections based on energy balance runs 
made during thermal investigations are also applied to each 
plate . 

Variations in the transpiration mass flux through each 
plate due to static pressure variations in the main air 
stream were found to be negligible in the present study. 

The largest pressure drop across the span of any plate, due 
to the main flow acceleration, is on the order of 10 percent 
of the pressure drop through the plate at the lowest blowing 
fraction of 0.001. For each static pressure distribution in 
the present experiments, thermal gradients, other than those 
observed by Moffat in zero pressure gradient flows, were 
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found not to exist in the plates when they are uniformly 
heated and blowing or suction is applied. 

D. Instrumentation 

The instrumentation used in the present experimental 
study is briefly described here. 

D. 1 . Temperature 

All measurements of gas temperature were made with 
calibrated iron-cons tantan thermocouples described in detail 
in reference [30], Temperature measurements relevant to this 
fluid dynamic study include free-stream static, porous plate, 
flowmeter, and ambient temperature. 

D.2. Injection and suction flow rates 

There are two flowmeters with different ranges 
for each of the twenty-four individual test plates. The 
meters, Fisher-Porter B4-27-10 and B5-27-10 Rotameters with 
float types BSVT-45-A and BSVT-64-A, respectively, provide 
measurement from 0.5 to 18.0 scfm per channel. Factory 
calibrations for these meters were checked against laminar 
flow elements, and against ASME orifice units. Agreement 
between these calibrations and the factory furnished calibra- 
tion curves was found to be within the accuracy of the test 
calibration; hence, the factory curves were used in all data 
reduction. Downstream of each flowmeter pair is a U-tube 
mercury manometer, measuring rotometer discharge pressure for 
the density correction to the indicated flow. 

D.3 . Pitot probes 

Three types of probes were used in the present 
investigation: flattened mouth pitot probes, Kiel probes, 

and Pitot- Prandtl probes. 

The flattened mouth pitot probes were used, for all 
boundary layer traverses. They are identical to those 


14 



utilized by Simpson [17] for tnis purpose and are described 
in detail by Simpson. Typical characteristics of tnese 
probes are 

1. Probe mouth dimensions; 0,012 incn by 0.035 inch 

2. Tube wall tnickness; 0.0025 inch 

3. Yaw and pitch plateaus: +10 degrees 

0.4. Static pressure ports 

Static pressure taps are located at 2 inch intervals 
along one side wall of the test section. They are spaced 
approximately 1 inch above the porous test surface so that 
they remain below the top cover for all test runs made in the 
present study, Tney are 0. 040-inch diameter square-edged 
holes flush witn the side wall with 1/8-inch O.D. brass tube 
hose connectors. According to Rotem [32], the static pressure 
error incurred by these taps, due to wall shear effect, is 
estimated to be about 0,0001 inches of water. 

These side wall taps were tested against the static 
pressure readings from Pitot- Prandtl probes located at the 
centerline of the test section. The static pressure ports 
of the probes were aligned with the side wall tap and then 
traversed through tne mainstream potential core. These tests 
were made with and without transpiration and for different 
upper wall settings to check for corner effects in the duct 
and streamline curvature in the potential core of the main 
air stream, respectively. The static pressure sensed by the 
Pitot-Prandtl probes was found to agree with that sensed by 
the side wall tap within 0.002 inches of water for all upper 
wall settings used in the present study. All present data 
were taken using the side wall pressure taps. 

D. 5 • Pressures 

Dynamic pressures were measured with four Dwyer 
model no. 100,5 inclined manometers with mirror scales. 
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reading directly in inches of water (-0.1 to 1.0) at 7b °F . 
Each manometer was periodically calibrated against a Harrison 
micromanometer, indicating a difference of no more than 
0.002 inches of water with 0,001 being the common value. 

D. 6. Distances 

Each probe was positioned relative co the porous 
wall with an individual traversing mechanism. These traversing 
mechanisms were mounted on a rigid support frame bolted to 
the structure supporting the test plate assembly. The longi- 
tudinal position of each traversing mechanism relative to a 
reference mark on the support frame is determined by means 
of spacers. Corrections were made for the longitudinal 
position of the probe mouth with respect to the traversing 
mechanism, yielding the distance from the beginning of the 
test section to the probe mouth. 

The micrometer-driven traverse mechanisms allow the probe 
to be advanced toward or away from the porous wall, with 
accurate alignment of the probe mouth in the flow direction. 
Although the least count is 0.001 inch, it was found that the 
micrometer could be advanced repeatedly to within 0.0002 inch, 
providing the micrometer was always advanced in the same 
direction to eliminate backlash. 

A more complete description of the traversing mechanism 
and the associated support frame is contained in reference 
[17]. 

E. Qualification of the Apparatus 

E. 1. Summary 

Velocity profile data taken on the apparatus were 
found to be "normal" and acceptable for incompressible flow 
over the impermeable flat plate for the range of free-stream 
velocities encountered in the present pressure gradient runs. 
Less than 3/8 percent variation was fov\ .i in the velocity in 
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the potential core. Discrepancies in experimental results 
using the two-dimensional energy integral equation were main- 
tained on the order of the experimental uncertainty of these 
data, suggesting chat if three-dimensional xlows exist in the 
present experiments they are weak. 

The test surface was found to he aerodynamically smooth 
for all impermeable flat plate data examined for the ran^e of 
free-stream velocities indicated. 

It appears that the test surface is aerodynamically smooth 
for all blowing, suction, and acceleration conditions investi- 
gated . 

E . 2 . Friction factors and mean velocity profiles for 
the impermeable flat plate 

Although pressure gradient flows are the major con- 
cern in the present study, the apparatus should produce accept- 
able impermeable flat plate data at the free-stream velocity 
levels encountered. The two requirements considered are: 

1. The generally accepted correlation of friction factors 

^ - 0.0128 Re 6 ~ 1/A (II-l) 

should be satisfied, as it was in the experiments of 
Simpson [17]. 

2. The profiles should exhibit wake strengths corresponding 
to those considered "normal" by Coles [33]. 

Impermeable flat plate runs were made at free-stream 
velocities of 42,86, and 126 ft/sec, respectively. An 
acceleration of the flow preceded tne region wnere the free- 
stream velocity was maintained constant. Comparison of the 
present data for U TO - 42 ft/sec with that of Simpson at the 
same level of free-stream velocity indicates that the effect 
of tr.'s initial acceleration is negligible. Figure 5 shows 
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the comparable agreement of these data with the friction 
factor correlation given by equation (II-l). The corresponding 
U + vs y + profiles are also shown to be in agreement with the 
accepted "law of the wall " in Figures 6 and 7 . 

Friction factors were obtained for these impermeable 
flat plate runs using a momentum integral equation method, 
described in Chapter III. Here, the concept of a "virtual" 
origin is utilized. The "virtual" origin is defined as that 
location from which the boundary layer appears to have begun 
assuming a uniform turbulent behavior everywhere. The 
"virtual" origin was calculated for each profile utilizing 
tne relation 

Y = 0.037 Re~° ’ 2 (II-2) 

This variation of 5p with X agrees with the friction factor 
correlation given by equation (II-l). The maximum percentage 
deviation of the "virtual" origin locations, corresponding to 
all profiles tor a given test run, is +3 percent. Friction 
factors were obtained by means of a power fit to the experi- 
mental Re 5 £ vs Re x variation, with X measured from the 
"virtual" origin location corresponding to the upstream pro- 
file for each test run. These friction factors are compared 
with equation (II-l) in Figure 5 > where the maximum deviation 
of the data from this relation is shown to be +7 percent. 
Inasmuch as the experimental uncertainty of chese data is 
approximately +5 percent, this agreement is considered accept- 
able for the purposes of the present study. 

Figures 6 through 9 show the mean velocity profile data 
in wall coordinates (U + vs y + ) . With the exception of the 
run lor = 126 ft/sec, the logarithmic region of the pro- 
files is best fitted by 

u+ = oTPT in |y + l + (II-3) 
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This uses the mixing length constant and intercept as proposed 
by Simpson. The result of using the values proposed by 
Coles [33] is also presented for purposes of comparison. The 
data for U^, = 126 ft/sec lie below both of these relations 
due to experimental uncertainties in C^/2 and possible sur- 
face roughness (see further discussion of surface roughness 
below) . 

Denoting the maximum deviation ~ from equation (II-3) 

U T 


as the Strength of the Wake , Coles [33] classified boundary 
layers as "normal" if their wake strength is within +20 percent 
of values given by the empirical relation 


r - 2 - 65 

T 


1-exp 


/-(Re 62 -500)\ 
[ 530 ) 


(n-4) 


for Re^ > 500 . This criteria results from Coles classifica- 
tion of nearly 500 impermeable flat plate profiles observed 
by independent investigators. 

The unblown flat plate data in the present study was 
tested for the Strength of the Wake . Using the log region 
fit given by equation (II-3), the Strength of the Wake values, 
given on Figures 6 through 9 > are within +10 percent of 
values predicted by equation (II-4). This agreement is 
similar to that observed by Simpson and, hence, the present 
impermeable flat plate data is considered "normal" and 
acceptable. 


E.3. Mainstream conditions 

Differential traverses made by Moffat [30] at axial 
locations along the duct show that, within the potential flow 
core, the velocity variation is less than 3/5 percent across 
any one plane. These measurements were made with the single 
screen system and no larger deviations were measured after 
the installation of the double screen system. 
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With the single screen system, hot wire anemometer 
measurements indicated a free-stream turbulent intensity 

vu*"2/U 00 of 1.2 percent for U^, = 44 ft/sec. The double 
screen system reduced this intensity level to 0.8 percent 
for Uoo = 44 and 25 ft/sec. A spectral analysis of this 
fluctuation data showed substantial contributions at the 
blower frequency (55 cps) and higher harmonics of this fre- 
quency. These may reflect the coupling of the blower noise 
to the tunnel system rather than a typically random turbulent 
fluctuation in velocity. Turbulence levels of 1 percent or 
more would normally result in distortions of the profiles 
from the "rormal" shapes described by Coles [33]. 

E. 4 . Two-dimensionality of boundary layers 

A unique method of determining the degree of two- 
dimensionality of the boundary layer is utilized in the 
present study. Direct measurements of Stanton numbers are 
available from which enthalpy thickness distributions can be 
calculated by integrating the two-dimensional energy equation.* 
These calculated values can be compared with those obtained 
using centerline velocity and temperature profiles at selec- 
tive stations. The discrepancy between the two estimates is 
indicative of the degree of two-dimensionality. 

For the most part, the calculated and measured enthalpy 
thicknesses agree within the experimental uncertainty of the 
data which is approximately +8 percent for all unblown and 
blown flows and which reach a maximum of +20 percent with 
strong suction (F = -0.002). Larger deviations are found 
to exist in some of the test runs for those profiles taken in 
the constant free-stream velocity approach region, but not in 
the pressure gradient region which is of major interest. This 


The reader is referred to reference [33] where all cor- 
responding thermal data is reported and discussed. 
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suggests that if weak three-dimensional flows exist in the 
present data they are acceptable in view of the experimental 
uncertainty of these data. 

For very low momentum and enthalpy thickness Reynolds 
number flows, small transverse variations in enthalpy thickness 
within the boundary layer, which are acceptable for thicker 
layers, become significant in this method of determining the 
degree of two-dimensionality. For this reason, it was found 
necessary to modify the apparatus for the test runs where 
K = 1.45 x 10" 6 in order to obtain the indicated agreement 
among estimates of enthalpy thickness. The installation of a 
double screen upstream of the nozzle and the addition of a 
trip inside the nozzle minimized the transverse variations 
which existed.* For an indication of the detailed character- 
istics of these variations, the reader is referred to 
reference [17] where transverse velocity and temperature profile 
data taken prior to this modification are presented. 

The discrepancies in the estimates of enthalpy thickness 
are given in Appendix A for each velocity profile where suf- 
ficient thermal data was available. With a few minor exceptions, 
this covers the velocity profile data corresponding to blowing 
and sucking fractions from F = -0.002 to 0.006. 

E.5. Surface conditions 

To satisfy the requirements of the ideal model, 

the porous surface must be aerodynamic ally smooth. In addition, 

the pore opening and pore spacing must be small such that the 

V inertia forces are small compared to the viscous forces at 
w 

the surface. 

Simpson [17] demonstrates, by means of appropriate models, 
that the latter of these two requirements is satisfied on the 
present apparatus. He further assumed that the plate "surface 


*A description of the screen and trip arrangements is given in 
Section B.l of this chapter. 


21 



*C -*«*y*‘i* ""T „£■«> **• ^ ... , ^W**"**. < * A*rt T *'•**’’* .'••■^^rF -*v1 


V 


x 



is in the plane of the particle crests for purposes of locating 
the "effective” no slip condition. The same location ^f the 
plate "surface" was assumed in the traverses of the present 
experimental study. 

Restricting attention to flows over an impermeable flat 
plate allows an estimate of the test surface roughness and 
its influence on the present velocity profile data. According 
to the data of Nikuradse [36], if the surface particles are 

, , . <lf k w 

within the viscous sublayer' ( ^ v — < 5 where k^ = particle 

size) the surface is aerodynamics lly smooth. Utilizing this 
criterion and the measured RMS roughness value of 0.00C2 

for k w , it follows, that for = 100 ft/sec and C^/2 = 

U k 

0.00300, ~ 0.6 and the surface is smooth. 

The impermeable flat plate data taken at different free- 
stream velocities yields a check on this preliminary analysis. 

In Figures 6, 7 , and 8, it is shown that for free-stream 
velocities of 42 and 86 ft/sec, U + vs y + relationships are 
the same in the logarithmic region. If roughness was important, 
it is expected that this U + vs y + relationship would change 
with . In Figure 9 > the profiles for U ro = 126 ft/sec 

are shown to be below the accepted "law of the wall" cor- 
responding to smooth wall data. This is thought to be the 
result of experimental uncertainties in C^/2 and possible 
surface roughness. 

It is evident from previous experiments [19-23] that 
acceleration has the effect of increasing the "viscous sub- 
layer" thickness of boundary layers on impermeable surfaces. 
Therefore, it is reasonable to conclude, on the basis of 
these results, that the test surface is aerodynamically smooth 
for all unblown experiments on this apparatus with and without 
acceleration where the level of is less than 100 ft/sec. 

With the exception of the impermeable flat plate run for 
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= 126 ft/sec, the level of is within this bound in 

all the present data. 

Simpson presents an argument concerning roughness effects 
with blowing and suction. Altnough the "viscous sublayer" 

is shown to decrease with blowing and increase with suction, 

U k 

the corresponding values of f° r strongest blowing 

and suction data are found to be small. For the highest 
suction case examined (C f /2 = 0.0076), it was found that 
U T k 

= 0.5 using k w = 0.0002 inch. For the highest blowing 

U k 

case examined, 1^/2 = 0.0001 and = 0.05 . These 

estimates, presented by Simpson, are based on U^, = 44 ft/sec, 
whereas the levels of found in the present experiments 

can be double this value. Taking this fact into account, the 

Ut kw 

corresponding estimates of yyy are still small enough to 

suggest that the test surface is aerodynamic ally smooth with 
blowing and suction. Sufficient proof of negligible rough- 
ness requires a direct basis of comparison which is not 
presently available in the cases of blowing and suction, and, 
therefore, this can only be taken as a plausibility argument 
as previously indicated by Simpson. 





Figure 1. Photograph of operating console area 


Figure 2. Photograph of test deck area 
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Figure 3. Test apparatus schematic 




Figure 4. Cross-section view of a typical compartment 
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Figure 5. C f /2 vs Reg 2 , unblown flat plate 




Figure 6. 

Unblown flat plate data of Simpson: run 3106 J, 

= 42 ft/sec 



Figure 7. 

Unblown flat plate: run 71967 , U^ = 42 ft/sec 
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CHAPTER II 


EXPERIMENTAL TECHNIQUES AND DATA REDUCTION METHODS 
A. Experimental Techniques 

The methods utilized to set up the desired flow conditions 
and the data taking procedures followed are described in this 
chapter. 

A . 1 . Setting of upper wall and transpiration rates 

The upper wall, the inlet velocity to the test 
section, and the transpiration rates must be carefully related 
to achieve an asymptotic boundary layer flow; i.e., one in 
which Peg 2 is constant in the flow direction. Ideally, the 
pressure gradient should be imposed at a position where the 
boundary layer already has the asymptotic value of Reg 2 • 

The mass flux at the wall must then be adjusted to maintain 
a constant blowing or sucking fraction F in the flow 
direction. These requirements were satisfied within reasonable 
ccuracy in the present experiments. 

The position where the pressure gradient should be im- 
posed was determined using available correlations for boundary 
layer flows in constant free-stream velocity. Simpson [17] 
has shown that, for unaccel rated ows where the mass flux 
at the wall is constant, 

c r d 

2 1 = c Re 5 ^ (HI-1) 


and 


H 


1 

1 “ (l + b) 1//2 + (l + o.635b ) 1/2 

to* ^ 


(III-2) 
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where c and d are unique functions of F . Applying 
these correlations to the asymptotic momentum integral equa- 
tion yields an equation which can he solved for the value of 
Re5 2 associated with given values of K and F , 

cRe 6 ^ « Res 2 (H + l) K - F . (III-3) 

Inasmuch as the boundary layer is allowed to develop with 
constant free-stream velocity in the initial portion of the 
test section, this yields an estimate of the position where 
the desired value of Res 2 will be attained. 

With the imposition of a strong pressure gradient the 
boundary layer is expected to adjust, resulting in different 
characteristics than those represented by the correlations 
utilized here. This, plus the fact that an abrupt acceleration 
is not possible, indicates that the above estimate is only 
approximate. For a given inlet free-stroam velocity the 
predicted position at which the flow is to be accelerated 
is quite insensitive to F , assuming that the "virtual 
origin" does not change for the unaccelerated portion of the 
boundary layer. Hence, acceleration was begun at the same 
position in the test section for all blowing and sucking 
fractions . 

As for the uppei wall setting in the pressure gradient 
region, a one -dimensional flow analysis indicates that an 
inclined flat upper wall, to a good approximation, yields a 
constant K flow for a given blowing or sucking fraction. 
Neglecting boundary layer displacement and streamline curvature, 
conservation of mass requires that the equation of the upper 
wall surface be given by 
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where 


KU (Xj 

= i + — [x-x a ] (m-5) 

Here, is the position where the pressure gradient is 

imposed. For all blowing and sucking fractions of concern 
this can be approximated by a linear variation within one 
percent. Based upon these considerations, an inclined flat 
upper wall was used to give a constant K flow. 

The specific top settings used are shown schematically 
in Figure 10. The inlet velocity to the test section and 
the slope of the upper wall in the pressure gradient region 
were chosen to maximize the length of the pressure gradient 
region and yet retain a turbulent boundary layer at the duct 
inlet . 

The transpiration rates needed to achieve a constant F 
for each test run were determined by means of an iterative 
process. The upper wall was first set for the desired com- 
bination of K and blowing or sucking fraction F . Tran- 
spiration flow rates were then sec up, based upon the free- 
stream velocity distribution of a zero transpiration test 
run. The transpiration flow rates were adjusted based on 
the resulting distribution of free-stream velocity. The 
process was continued until the desired blowing fraction had 
been attained within the desired limits. 

A. 2. General test procedure 

A specific test procedure was followed during each 
test run. After the air system had stabilized, the flow- 
meters were given final adjustments and the thermal data 
were taken. The traversing instruments were then positioned 
using spacers and. a reference mark on the traversing support 
beam. Free-stream pitot readings for each boundary layer 
probe and a Kiel probe, all placed in the constant free-stream 
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Figure 10. 

|| Duct conf igurations for pressure gradient runs 
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velocity region of the channel, were compared to check for 
line leaks and possible fouling of the mouth of the boundary 
layer probes. After these initial adjustments and checks 
were made, velocity traverses were taken according to the 
procedure described below. When approximately half of the 
traverses were obtained, the static pressure distribution 
along the channel was measured and recorded along with the 
flowmeter settings and thermal data. All test conditions 
were monitored at intervals throughout the entire test run. 

A. 3 . Traversing procedure 

Velocity traverses were taken on the duct centerline 
in a manner similar to that used by Simpson [17]. The boundary 
layer probe was advanced toward the wall until contact be- 
tween the mouth of the probe and the wall was visually evident. 
The probe was then backed away from the wall in 0.001 inch 
increments according to the micrometer scale. More than one 
reading was obtained with the probe on the wall during this 
procedure. The position where the probe left the wall was 
determined by means of a dynamic pressure method. This method 
is based on the fact that the indicated dynamic pressure re- 
mains the same as long as the probe touches the wall but in- 
creases markedly as it leaves. According to Simpson, this 
method of wall location agrees with an electrical contact 
method within 0.0005 inches. The probe was then advanced 
through the boundary layer in approximately °qual intervals 
of dynamic pressure. 

All dynamic pressure readings were obtained with the 
instrumentation described in Chapter II. The side static 
pressure ports were used as references in regions of the 
duct where the static pressure exceeded the 1-inch of water 
limit of the Dwyer manometers. The local static pressure 
necessary for the calculation of velocities was determined 
by linearly interpolating between static pressure ports 
adjacent to the mouth of the probe. The use of side ports 
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as references has the advantage of reducing the effects of 
minor time-variant fluctuations due to blower and ambient 
changes. The response of the probe-manometer system to 
within 0.001 inch of water of the approached value is estimated 
to be approximately 30 seconds but readings were taken over 
an additional one-minute interval to obtain a time average. 

Prior to each traverse, manometers were zeroed and 
test conditions were monitored. These checks were also 
made at the completion of each traverse as a check on system 
changes during the data taking period. 

B. Data Reduction 

The raw experimental data was initially reduced by means 
of a computer program. A statement of the program in Fortran 
IV language is given in Appendix B. In its most general terms, 
the program performs the following functions: 

1. Accepts input and converts to standard units. 

2. Corrects for calibration curves of instrumentation. 

3. Calculates velocity profile integral parameters using 
the trapezoidal rule. 

4. Calculates dimensionless parameters. 

5. Summarizes and prints out a raw data listing and 

* 

a summary of output data with uncertainty intervals. 

In addition to the above, the output data were punched on 
computer cards which were utilized as input cards to secondary 
computer programs. The secondary programs were either cor- 
relation reduction programs or designed to produce plots and 
the final data tables, presented in Appendix A. 

Specific features of the data reduction process are of 
considerable importance in the interpretation and correlation 


The uncertainty analysis is not presented in the program 
listing for purposes of clarity. 
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of results. For this reason, they are described in detail 
here. The method of reducing pitot tube readings and the 
associated uncertainty analysis are discussed. The method 
utilized to determine local blowing rates and pressure 
gradients is presented. The means utilized to obtain experi- 
mental estimates of friction factors are also presented in 
detail. 


B . 1 . The reduction of pitot tube readings 

The pitot tube data were reduced to values of mean 
velocity by application of the Bernoulli relation 


U 



2 ( p „- p „)s r 


p 


(III-6) 


after suitable correct ons were applied. The corrections are 
applied, to account for errors inherent to the flattened -mouth 
probes utilized. 

Simpson [17] considered the importance of turbulent 
fluctuation, viscous, shear, yaw and pitch, and wall effects 
on pitot tubes of this type. He concluded that: 


1. No correction for yaw or pitch angle need be applied 
for all the blowing and sucking cases he studied. 

He found that the mean streamline angles, tan -1 (v w /u)> 
were less than 10 degrees which is within the pitch 
plateau experimentally determined, for these probes. 

Yaw angle corrections could be avoided by positioning 
the probe accurately. 


2. 


The viscous correction that should be applied, is 
based upon the data of Hurd., Chesky and Shapiro [34 ] . 
Letting 



2( p s - p 0 )g 


c 
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the correction he proposed, is given by 


o » 1.02 - — 5 — + 2^51 (m-8) 

Re p 

for I\t,p < 40 and C p = l.C for Re p > 40 . 

The Reynolds number Re ' is based on the outside 
hydraulic diameter of the probe mouth and the local 
velocity. 

3. The turbulent fluctuation components of velocity 
can be neglected, resulting in only slight errors 
in the uncorrected velocity profile. 

4. Negligible shear gradient and wall effects were 
found in the severe case of an asymptotic suction 
layer. Therefore, these effects were not considered 
in the remaining cases he studied . 

The same set of corrections is used in the present study. 
The mean streamline angles are found to be within the bounds 
stated by Simpson and, hence, yaw and pitch angle corrections 
are not applied. Since turbulent fluctuation components were 
not measured, a turbulent fluctuation correction cannot be 
applier" . Shear gradient and wall effects are neglected 
inasmuch as the cases studied here are not as severe as the 
asymptotic suction layer flow studied, by Simpson. 

B.2. Determination of local blowing rates and pressure 
gradients 

As indicated in Chapter II, the mass flux through 
each plate can be considered uniform within the limits of the 
plate porosity variation. The effects of the static pressure 
drops encountered in the mainstream flow direction for those 
cases considered here were found to be negligible. Hence, 
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the mass flux at a given velocity traverse station is taken 
to be that corresponding to the center 6-inch span of the 
associated plate. The blowing fraction F , in turn, is 
based upon this mass flux and the local free-stream velocity 
for each traverse. The same approach is taken for the suction 
cases. 

The fact that constant K flows are being studied is 
utilized to determine local values of the applied pressure 
gradient. The pressure gradient at each static port is 
taken to be that corresponding to a linear velocity distribu- 
tion between the adjacent static ports. Selected values of 
thf» resulting pressure gradients are found to agree with 
values obtained, from a plot of the free-stream velocity dis- 
tributions, within +5 percent. The free-stream velocities 
at each static port station are determined from the free-stream 
velocity head, measured in the constant free-stream velocity 
approach region of the duct, and the measured static pressure 
distribution. These measurements are sufficient for this 
purpose only if the total pressure head remains constant and, 
hence, a potential core was maintained, during each test run. 
These known quantities at the static port stations allow local 
pressure gradients to be determined by linear interpolation 
of the parameter K between adjacent static ports. 

B. 3 . Experimental determination of friction factor 

The friction factor C^/2 is experimentally deter- 
mined. by four different methods. These methods are of the 
indirect type and all four are not applied to the same profile. 
Each of these methods is described, with assumptions stated 
and an indication of where each is applied in the present 
study. Due to the propagation of experimental errors in the 
use of these methods, it is found necessary to obtain a "best 
estimate" of friction factor. The approach taken to obtain 
this estimate is also described. 
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B.3.a Momentum integral methods 


Two methods of obtaining friction factors 
from the two-dimensional, constant property momentum integral 
equation are utilized. One method is applicable in boundary 
layer flows where the free-stream velocity is maintained 
constant, whereas the other is only applicable in near- 
asymptotic boundary layer flows of the type considered in 
the present study. 


i . Constant free-stream velocity 

Restricting attention to two-dimensional 
boundary layer flows over impermeable flat plates, the constant 
property momentum integral equation can be expressed 

C f dRegp 

r- = 3 rS^ (hi-9) 

X 

With measured values of momentum thickness at successive 
X-stations, this relation can be used to obtain values of 
Cf/2 for such flows. In order to reduce the error inherent 
in the evaluation of the derivative involved, the following 
power fit is assumed to apply. 


Re 62 = aRe* (IU-10) 

For Re measured from the "virtual origin", this is considered 
representative of the actual flow condition. Substitution of 
this relation into the above momentum integral equation yields 
the following expression for the friction factoi C^/2 , 

C- b-1 

= abRe x (III-ll) 
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Hence, estimates of C^/2 are obtained from an experimental 
determination ol the constants a and b . A similar approach 
can be taken in cases where either blowing or sucking at the 
wall exists. The reader is referred to reference [17] for a 
discussion of this more general approach. 

This particular method requires more than one velocity 
traverse as well as a determination of a "virtual origin". 

For this reason, it was only utilized for the surface rough- 
ness investigation discussed in Chapter II where both require- 
ments were satisfied. 

ii. Near-asymptotic flow 

As indicated in Chapter I, the constant 
property two-dimensional momentum integral equation can be 
expressed as 


C r dRe&p 

= Reg 2 (H + 1)K - F + ■ ■ . (III-12) 

^ x 

For near-asymptotic flows such as those considered in the 
present study, 

dRe&p 

- < < Re 62 (H + 1)K - F (III-?3) 

X 

This characteristic allows equati >n (III-12) to be used to 
calculate friction factors in this type of flow on the basis 
of local experimental data with a small correction term, 
d.Re5 2 

-rp , applied. The relative magnitude of the correction 

term does increase with an increase in blowing fraction. This 
is the result of the right side of equation (III-13) tending 
toward zero with an increase in blowing fraction. Similar 
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undesirable effects are present in all momentum integral 

approaches when blowing is considered. 

The following procedure is followed in the evaluation of 

dRe6 2 

the correction term . In terms of the local X-Reynolds 

X 

number, this term may be expressed 

dRe & 

’ BRi-^ 15 Re x + l) 011-141 

X 

This relation is used with an assumed fit to the experimental 
data given by 



Reg 2 ~ m exp[n Re x ] (III-15) 

where the experimentally, c-termined constants m and n are 
based on successive velocity profiles taken in the pressure 
gradient region of the duct. If possible, the profile down- 
stream of the X* station of interest was used for this purpose. 
The term was evaluated on the basis of the local X-Reynolds 
number rather than the integrated X-Reynolds number to avoid 
a propagation of uncertainties in the pressure gradient 
parameter K . Equation (III-15) is found to fit the trend 
of the c.xporimental data. 

Equation (III-12) was used to obtain values of C^/2 
for all profiles measured in the pressure gradient region of 
the duct. For purposes of confirmation and slight adjust- 
ments, the viscous sublayer model method described in the 
following section was used as an independent means of obtain- 
ing values of c^/2 . 

B.3.b Viscous sublayer model method 

This method relies on tjio fact that in a 
thin region necr the wall uolecular viscosity' governs the 
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flow. Neglecting X-derivatives and treating the flow as 
laminar in this region, the constant property boundary layer 
equations yield the solutions 


U_ 

U„ 

oc 


1 

F 





for F ^ 0.0 


(III-16) 


and 



for F = 0.0 . (III-17) 


These relations can be solved for C f /2 , allowing the cal- 
culation of friction factors from experimental values of 
V , U , , and y . Here, U must correspond to a y- 

position within the region of application of these solutions. 

The reliability of this method is dependent upon the 
ability to measure low dynamic pressures within the "viscous 
sublayer" and the accuracy within which the probe position 
relative to the wall is known. Hence, it was used to confirm 
friction factors obtained by means of the momentum integral 
equation. Friction factors corresponding to profiles 
traversed at stations in the entrance and recovery sections 
of the duct were also estimated by this method. 

B.3.c Logarithmic region method 

This method of determining friction factors 
> is based on the assumption that a "universal law of the wall" 
exists. The specific form of this law proposed by Simpson 
[17] for the impermeable wall case is given by 
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07 W Jn|y + | + 5.55 


(III- 18) 


U 


+ 


This relation is assumed to apply to those regions of the duct 
where free-stream velocity is maintained constant for all un- 
blown test runs. Velocity profile data can be cross-plotted 
against this assumed relation to yield an estimate of the 
associated friction factor. 

Various "laws of the wall" for turbulent boundary layers 
have been proposed [15>17>29] which account for the effect 
of blowing and applied pressure gradients. Since these are 
in disagreement with one another and are not based on data 
for the combined problem, they are not considered in the 
present evaluation of friction factors. 

B.3.d Determination of best estimate of friction 
factor 

The "best estimates" of friction factor were 
obtained by blending the experimental velocity profiles with 
the viscous sublayer equation. 

The graphical procedure used in the pressure gradient 
region was based on the following observations: 

1. The viscous sublayer equation was not sensitive to 
small changes in the friction factor so that the 
momentum integral values could be used to predict 
sublayer profiles. 

2. On the basis of the viscous sublayer equation and 
the U + vs y + profiles corresponding to the momentum 
integral values of C^/2 , the inner regions of 

the sublayer (y + < 15) appeared invariant in the 
streamwise direction for the majority of the data. 

For a given test run, all profiles in the pressure gradient 
region were inspected simultaneously. The profile which 
best matched the viscous sublayer equation was selected as a 
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reference profile. A slight adjustment in the value of 
Cj/2 for the reference profile was made for some of the 
test runs if agreement with the viscous sublayer equation 
was not obtained to the degree found in the majority of the 
data. The remaining profiles were then matched to the viscous 
sublayer equation with the requirement that the velocity 
distributions in the inner regions (y + < 15) correspond to 
that for the reference profile. 

In this procedure, a first estimate of friction factor 
was obtained by the momentum integral method. The resulting 
"best estimates" of C f /2 were found to agree with the 
momentum integral values within the uncertainty of the cor- 


rection term 


dReB 2 

dRe~ 


(Re x K 


applied in the momentum integral 


method. The calculated uncertainties in the correction term 
represent variations in friction factor of approximately 
+0.00020 , +0.00030 and +0.00040 in the present data for 
K = 0.57 x 10"6 t 0.77 x 10 “6 and 1.45 x 10~6 , respectively. 
This amounts to approximately 15 percent of friction factor 
for the unblown runs and to as much as 50 percent for the 
highly blown runs. 

For the unblown test runs, the logarithmic region method 
values of friction factor were considered "best estimates" for 
all constant free-stream velocity profiles. 

For the blowing or suction cases, "best estimates" of 
friction factor were found for all constant free-stream 
velocity profiles by graphically matching the U + vs y + 
profiles to the viscous sublayer equation. If data points 
were not found to exist in the viscous sublayer region an 
extrapolation of the existing data points was used. This 
extrapolation was based on experimental profiles of Simpson 
[17] corresponding to approximately the same values of Re§ 2 
and F . In this graphical matching procedure as many data 
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points as possible were used and not just two data points 
near the wall as was done in the viscous sublayer method. 

The deviations of the other experimental determinations 
from these "best estimates" are discussed in Chapter IV with 
other experimental results. 

B. 4. Uncertainty intervals 

Only those uncertainties resulting from the necessary 
interpolation between divisions on an instrument and those due 
to fluctuating values of the measurand are considered in the 
present uncertainty analysis. Errors in instrumentation 
calibration are not considered since they are "fixed errors" 
and uncertainties due to uncontrolled peripheral variables 
of the experiment cannot adequately be estimated. 

The procedure of Kline and McClintock [37] is followed 
to account for the propagation of uncertainties in all cal- 
culated quantities. The following basic uncertainty intervals 
are assumed with 20:1 odds: 

Distance from the wall 

Flowmeter reading (25 cm 
full scale) 

Dynamic pressure 

Relative locations of 
static ports 

Location of probe relative 
to static ports 

The resulting estimates of uncertainties in the calculated, 
quantities are indicated in the data tabulations found in 
Appendix A. 


0.0005 inches 
0.02 centimeters 

0.002 inches of water 
0.008 inches 

0.0l6 inches 
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C. Calculation of Shear Stress Profiles 


Shear stress profiles were generated from the measured 
mean velocity profiles in order to check commonly assumed 
shear stress distributions in the inner regions of the boundary 
layer. The relation used for this purpose is developed here. 

The following assumptions are made: 

(1) U/U^ vs y/6 similarity holds for the entire 
boundary layer, i.e., X dependency contained in 
6 , where 6/6g = constant. 

(2) The contribution due to the non-asymptotic condition 
dReSo 

-^5 ^ 0 is small in the inner regions. 

aR x 

The importance and applicability of these assumptions is 
presented in the following development and discussion. 

For constant property two-dimensional boundary layers, 
the X-momentum and continuity equations can be expressed in 
the integral forms 


and 



d_ 

dx 


rV 


U dy + UV 


(III-19) 


V 



d_ 

dx 


y 

U dy , 


o 


(III-20) 


respectively. Substituting equation (III-20) into equation 
(III-19) yields 
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u dy 

(III-21) 


By means of simple operations of calculus and appropriate 
definitions, the following identities are obtained. 


d 

dx 


U 2 dy 


6U 2 5- 
00 dx 


.y/6. 2 

(i) d(y/6) 


db 0 dU ~ 
5 dx + 2 U ~ dx 


•y 


fej 


dy 


(III-22) 


and 


k \ u d y - 6 «. 


.y/B 


00 dx 


Gfc) d(y/6) 


IJ dU m 

oo q 5 oo 

b - dx + dx 


(»J 


dy 


(III-23) 


Substituting these identities into the right-hand side of 
equation (III-21) and. applying the similarity assumption (l) 
gives 
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00 do TI “ 

6 dx 00 dx 
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(III-24) 


For 5/82 = constant , a result of assumption (l), it can be 
shown that 



d6 

dx 


U^ 2 dRe 62 dU^ 

6 2 ^ R x ” U °° 


(111-25) 


Substitution of equation (III-25) into equation (III-24) 
yields 



(III-26) 

Dividing both sides of equation (III-26) by t w , transposing 

+ -f + 4* 

terms, and applying the definitions of t , U , y , V w , 
p + , and Cf , the following expression is obtained. 
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(111-27) 

This was used for calculating the distribution of shear 
stress. The integrals involved were evaluated using the 

dRegp 

trapezoidal rule and the derivative -rs was taken to be 

dR x 

that used in the momentum integral method of estimating 
friction factor. 

Shear stress profiles were computed for traverses in a 
region of the duct where profile similarity was found to 
exist in accordance witn assumption (l) . Assumption (2) 
does not explicitly appear in this development but limited 
errors due to uncertainties involved in the evaluation of 

dRegp 

the derivative ^ . 
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CHAPTER IV 

THE EXPERIMENTAL RESULTS 

The experimental data consist of mean velocity profiles 
obtained, in near-asymptotic boundary layer flows where the 
pressure gradient parameter K and blowing fraction F are 
maintained constant. Three pressure gradients were invest!- 
gated: K = 0.57 x 10 , 0.77 x 10 , and. 1.45 x 10" D . For 

each pressure gradient, the conditions investigated, cover a 
range of . niform blowing fractions frc-’' F = -0.004 to 0.006 . 

The estimates of friction factor and. their consistency 
and reliability are discussed in this chapter. Mean velocity 
profiles are also presented and. checked, for profile similarity 
and development characteristics. 

The asymptotic characteristics of these boundary layer 
flows are discussed. Shear stress distributions, computed 
for purposes of representing sublayer data, are also presented. 

All data presented graphically are also tabulated, in 
Appendix A with estimates of their uncertainties. 

A. Friction Factor Data 

Estimates of friction factors were obtained, for each of 
the experimental velocity profiles. The methods utilized in 
obtaining these estimates are discussed, fully in Chapter II 
and. the friction factor data is presented and discussed here. 

Restricting attention to profiles taken in the pressure 
gradient region of the duct, two experimental determinations 
of C f /2 are available. In Appendix A, the values of 
Cf/2 given correspond to: 

1. Near-asymptotic momentum integral method 

2. Viscous sublayer model method (average of values 
corresponding to first two profile data points) 

3. "Best estimates" of C^/2 utilized in succeeding 
data reduction. 
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For 63 out of a total of 79 profiles) the ..amentum integral 
values of friction factor agree with the "best estimates" 
within +10 percent. For 6l profiles, the viscous sublayer 
values are also found to agree with the "best estimates" 
within +10 percent. In the majority of these profiles where 
agreement was not found with the viscous sublayer method, 
the "viscous sublayer" had not been adequately penetrated 
during the traversing procedure. Only 2 profiles are presented, 
where the momentum integral value and viscous sublayer value 
do not agree with the "best estimate" within this +10 percent 
range . * 

In the constant free-stream velocity approach and recovery 
regions of the duct, only the viscous sublayer method was used. 
One exception is the impermeable wall case where the logarithmic 
region method was used. Slightly less than half of the profiles 
in this region had agreement between the viscous sublayer 
values and "best estimates" of C f /2 within +10 percent. This 
poor agreement is explained, for the most part, by the inability 
to penetrate the "viscous sublayer", a predominant characteristic 
in the recovery region of the duct. 

Simpson [17] correlated friction factor data on the basis 

cf the blowing parameter B = Q-jg and the unblown friction 
factor correlation 


Uf o -I/ 2 * 

— = 0.0130 Re t2 . (IV-1) 

He found that, for constant free-stream velocity and a 
variety of blowing and succion distributions, the friction 
factor data obtained on the present apparatus could be fitted 


These profiles are for run 81668 , M * 4, and 5 . 
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by the relation 



(IV- 2) 


The friction factor data of Simpson, for constant and slowly 
varying A" cases, agreed with equation (IV-2) within his 
calculated experimental uncertainty. Simpson sets this un- 
certainty at approximately +10 percent for B <0 and as 
much as + 25 percent for the higher values of B encountered 
in the present study. 

The "best estimates" of friction factors for the present 
pressure gradient data and constant free-stream velocity data 
are compared with equation (IV-2) in Figure 11. The constant 
free-stream velocity data a 3 shown to be consistent with this 
correlation within the calculated uncertainty and no anomalous 
effects are noted. Friction factors for the pressure gradient 
data are shown to lie 10 to 20 percent above the values given 
by equation (IV-2) for all accelerations combined with blowing 
and below by approximately 10 to 25 percent for all accelera- 
tions in the case of suction. 

It is concluded that the "best estimates" of friction 
factor for the present velocity profile data are self-consistent 
and are in acceptable agreement with the experimental deter- 
mination according to either the momentum integral method or 
the viscous sublayer method. It is found that the present 
friction factor data for constant free-stream velocity agrees 
with that of Simpson. The effect of acceleration on C^/2 
is found to differ in the cases of blowing and suction. In 
either case this results in relatively small deviations from 
equation (IV-2). 
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B. Mean Velocity Profile Data 


The characteristics of the mean velocity profile data 
are discussed! here. Profile similarity in both the inner 
regions and outer regions of the boundary layer are considered. 
The development characteristics and asymptotic behavior of the 
boundary layer flows are also presented. 

B.l. Inner region development and similarity 

Velocity profiles are present in wall coordinates 
(U + vs y + ) in Figures 12 through 20 utilizing the "best 
estimates" of friction factor. For purposes of comparison, 
the accepted "law of the wall" with constants proposed by 
Simpson [17] is also presented on each of the graphs. The 
profile obtained in the constant free-stream velocity approach 
region is presented along with the profiles obtained in the 
pressure gradient region of the duct. 

For the impermeable wall case, all three pressure 
gradient runs are presented. It is shown in Figures 12 
through 14 that the inner region of the boundary layer respond 
rapidly to the imposed pressure gradient and assume a unique 
distribution corresponding to a given value of K . For the 
strongest pressure gradient, K = 1.45 x 10 - ^ , slight adjust- 
ments of the layer are found to exist through the entire 
pressure gradient region. It is concluded that, in the im- 
permeable wall case, similar inner region profiles exist in 
these near-asymptotic boundary layer flows, and. that the 
shape of the profile is dependent upon the value of the local 
pressure gradient parameter K . 

Two characteristics of these boundary layer flows are 
shown in these inner region coordinates: (l) The profiles 

deviate from the acceptea "law of the wall" in an "overshot" 
manner within the logarithmic region and (2) the wake region 
is substantially diminished. The indicated profile "over- 
shoot" is specifically an upward, displacement of the U 4- vs y + 
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profile from the "flat plate law of the wall" in the fully 
turbulent region of the layer. Similar qualitative charact- 
eristics have been observed by Launder and Stinchombe [20], 
and Patel and Head [24], Patel and Head reported an initial 
"undershoot" of the logarithmic region, not seen in the 
present data. In the present study, traverses were not taken 
at stations along the duct where initial adjustments of the 
layer to the imposed pressure gradient occurred. Hence, the 
indicated results do not exclude this possibility. For the 
present data, it is shown that the degree of "overshoot" in 
the logarithmic region increases with K . This behavior can 
be represented as an increase in the thickness of the "viscous 
sublayer" region. The diminished wake is a direct result of 
the low shear stress in the outer regions of the layer, a 
characteristic associated with favorable pressure gradients. 
The resulting effect of these structural changes is to produce 
an "overshoot" profile in the logarithmic region where the 
logarithmic region decreases in its extent with the parameter 
K . This behavior is in agreement with the continuous re- 
laminar izat ion observed by Patel and Head in accelerated 
flows stronger than those in the present study. 

The relative constancy of the friction factors for these 
unblown constant K flows is also indicated in Figures 12 
through 14. In inner region coordinates the natural pressure 


gradient parameter is 


p + = 


Sc*S. 


dP 


according to a dimen- 


? dx * 

'w "1 

sionless momentum equation. For constant property flows, 
this parameter can be expressed as p + =-K/(Cf/2)3/2 # Hence, 
these constant K boundary layer flows are also characteristic 
of constant p + flows. The inner regions are, therefore, 
related on the basis of this parameter for use in the theoret- 
ical prediction of results. The theoretical predictions are 
presented in Chapter V. 
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In Figures 15 through 18 , similar effects of acceleration 
are shown to exist with blowing at the wall. The adjustment 
of the layer in the inner regions is found to be greater from 
station to staoion with blowing. The degree of overshoot the 
profiles exhibit above the flat plate profile is found to be 
decreased with blowing. This behavior is possibly the result 
of the increased turbulence level in the boundary layer due to 
blowing. The wake region is also found to be decreased to a 
greater extent than in the unblown layer, indicating a greater 
increase in friction factor with acceleration. 

These similar inner regions of the boundary layer are 
related on the basis of the natural parameters p + and V w + 
in Chapter V. 

The boundary layer flows in the case of suction at the 
wall are shown in Figures 19 and 20 to have similar charact- 
eristics as found in the unblown and blown layers. The profile 
"overshoot" is found to be greater as expected. For moderate 
suction the adjustment of the layer to one exhibiting similar 
inner regions is shown to exist except in the case of a very 
strong acceleration, where the structure of the inner region 
attains predominately laminar characteristics (roundness of 
profile). In all combined pressure gradient and suction runs, 
substantial structural changes of this type in the inner region 
were noted with the exception of those flows where 
K = 0.55 x 10“ 6 , and F * - 0.001 and - 0 . 002 . 

It is concluded that the qualitative characteristics of 
boundary layer flows for acct leration apply for the range of 
blowing and sucking fraction;? considered in the present study. 
The quantitative effects are dependent on blowing or sucking 
fraction and similar-profiles in the inner regions are found 
to exist except in flows where the structure of the inner 
regions of the layer appears to be substantially and continually 
changing, which suggests that relaminarization is perhaps 
occurring. 
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B.2. Outer region development and similarity 

The development of the outer regions of these 
asymptotic and near-asymptotic boundary layer flows is shown 
in Figures 20 through 26. Here U/U,,,, vs y/6 profiles are 
presented comparing the profile in the approach region with 
those obtained in the pressure gradient region of the duct. 

It is demonstrated that, for all blowing and sucking 
fractions considered, the layers approach and attain an 
asymptotic limit, characterized by similar profiles in these 
coordinates. The flows where the sucking fraction was 
F = -0.004 are exceptions to this rule since the profiles 
are found to continually adjust to a laminar mode, a char- 
acteristic not indicated in the flows graphically represented. 
It is noteworthy that similarity in these coordinates also 
implies similarity in "velocity-defect*' coordinates since 
friction factors in the flow direction are found to be 
essentially constant.* 

In the impermeable wall case, this similarity behavior 
indicates that asymptotic boundary layer flows for the present 
range of K are also constant P flows where P < -0.5 , 
as shown in Chapter I. 

The response of the layer to its asymptotic condition is 
shown to be extended farther down the pressure gradient region 
of the duct with an increase in blowing fraction, which is 
reasonable in view of the greater adjustment found necessary 
with blowing. The suction run with the strong pressure 
gradient, K = 1.45 x 10”^ , is shown in Figure 27 not to 
follow this trend, suggesting structural changes in the 
boundary layer as noted in the discussion of the inner regions. 

It is concluded that in the asymptotic and near-asymptotic 
boundary layer flows investigated in the present study similar 


* "Velocity-defect" coordinates are -w — vs y/6 . 
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profiles were attained in the outer regions of the layer, 
except in the flows where the sucking fraction was -0.004 , 

This coupled with a similar conclusion relating to the inner 
region of the layer confirms the existence of completely 
similar profiles in asymptotic boundary layer flows, with the 
exception of the suction cases noted in the discussion of 
the inner regions. 

B.3. Asymptotic values of local momentum thickness 
Reynolds number and shape factor 

By definition the asymptotic boundary layer flows 
considered in the present study are characterized by a constant 
value of Re6 2 • The present experimental data indicates that 
unique values of Regg and shape factor H correspond to 
given values of the pressure gradient parameter K and 
blowing fraction F . 

In Figure 28, the values of Re§ 2 » corresponding to the 
last two or three profiles obtained in the pressure gradient 
region of the duct are plotted against K , for parametric 
values of F . The agreement between the values of Reg 2 
associated with consecutive profiles for a given run shows 
the success in attaining an asymptotic condition for all 
three values of K where F > -0.002 . The continuous trend 
of the data displayed in this figure indicates these asymptotic 
conditions are unique in view of the fact that the values of 
K correspond to different ranges of free-stream velocity in 
the present experiments. 

In Figure 29, the experimental estimates of the asymptotic 
values of Re6 2 for the impermeable wall flows are compared 
with those values corresponding to an assumed turbulent 
velocity profile and an exact laminar solution. The assumed 
turbulent velocity profile is given by the equation 

U + » 8.7(y + ) 1/7 (IV-3) 
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Applying this equation to the asymptotic form of the two- 
dimensional momentum integral equation yields the relation 


Re 52 | 

|asymp 


0.0160 K“ 0,8 . 


(IV-4) 


An exact solution is available for the case of laminar flow 
[28], Applying this to the asymptotic form of the two- 
dimensional momentum integral equation yields the relation 



asymp 


0.3467 K' 0,5 


(IV-5) 


It is shown that the values predicted by the 1/7-power profile 
assumption are in agreement with the experimental values, 
whereas, the laminar values are considerably lower. On the 
basis of these results, unblown asymptotic boundary layer 
flows for values of K as high as 1.45 x 10~ 8 are believed 
to be turbulent in nature. 

The shape factor corresponding to profiles in the pressure 
gradient region for a given run are found to attain a reasonably 
constant asymptotic value with the exception of the strong 
suction cases for F = -0.004 . 

In an attempt to obtain a correlation of shape factor 
with pressure gradient, the experimentally determined values 
are compared in Figure 30 with the shape factor correlation 
proposed by Simpson [17] for turbulent boundary layers. This 
correlation of constant free-stream velocity data is given by 


H 



(IV- 6) 
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where 


C = 



CV2 


(1 + B) 1/2 + (1 + 0.635B) 1/2 


I 


(IV- 7) 


The present constant free-stream velocity data is also shown 
on this figure, indicating that this equation represents a 
good mean of the zero pressure gradient data even though the 
data scatter is relatively large. 

No significant correlation of shape factor with pressure 
gradient can be determined. The experimental shape factors 
also appear to be reasonably insensitive to blowing fraction 
for values of F greater than or equal to -0.001. In the 
cases of stronger suction, the values of shape factor are 
found to be higher than those indicated in Figure 30 and are 
not shown due to the strong possibility that these boundary 
layers do not behave in a turbulent manner. It is also note- 
worthy that the shape factor correlation of Simpson is shown 
to yield higher values of shape , factor than the experimental 
values presented and is, therefore, not considered applicable 
in such strong favorable pressure gradient flows. 

The shape factor corresponding to the exact laminar 
solution for the unblown asymptotic layer has a constant 
value of 2.0 . In the present unblown data, values of shape 
factor on the order of 1.3 correspond to those profiles in 
the pressure gradient region of the duct. Hence, the unblown 
asymptotic boundary layer flows considered are believed to be 
turbulent in nature, on the basis of shape factor as well as 

Re 62 • 

A similar comparison of Resg and shape factor H is 
not presented for the blown and sucked flows since exact 
laminar solutions do not exist in these cases. However, it 
is reasonable to conclude that the present blown boundary 
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layer flows are turbulent in nature, although transition to 
a laminar mode may occur for strong suction. 

C. Shear Stress Profiles 

In formulating models that represent the present mean 
velocity profile data, it is desirable to know the distribu- 
tion of shear stress through the boundary layer. This knowledge 
may allow representation of the data using simple eddy-viscosity 
models. Neglecting X-derivatives in the inner regions of the 
layer, the X-direction momentum equation yields the shear stress 
distribution given by 


t + = 1 + u + V * + p + y + • (rv-8) 

w 


Here, the dimensionless shear stress t + is defined as 

t + « — . This distribution has been proposed by numerous 

T w 

investigators as representative of that in the inner regions 
for small pressure gradients. In the present study shear 
stress profiles were computed to test the applicability of 
equation (IV-8) in strong pressure gradients. 

The following relation, developed in Chapter III, 
allowed the generation of shear stress profiles from the 
present mean velocity data where similarity in U/U^ vs y/6 
coordinates existed. 


+ 


1 + u + v w + + y + p + 



1 

c f 

62 2 1 


dRe6 2 

d^T" 





+ 



(XV-9) 
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The resulting shear stress profiles are tabulated in Appendix 
A. 

In Figures 31 and 32, representative shear stress profiles 
are compared to those obtained utilizing equation (IV-8) for 
the pressure gradients K * 0.57 x 10"^ and K » 1,45 x 10"^ , 
respectively. It is observed that the maximum shear stress 
occurs, in all profiles presented, for values of y + less 
than 25. It is also shown that, beyond the maximum shear 
position, the experimentally determined profiles substantially 
and increasingly deviate from that given by equation (IV-8), 
Inasmuch as the inner regions of the layer exhibit this 
behavior, the shear distribution given by equation (IV-8) 
is found not to be applicable in semi-empirical representations 
of the experimentally determined profiles. 

Upon comparing the relative contributions of the terms 
in equation (IV-9)> it is found that up to y + = 140 the 
contribution of the non-asymptotic term is less than 2 percent. 
This results from the fact that asymptotic or near-asymptotic 
flows are considered in the present study. Neglecting the 
non-asymptotic term in equation (IV-9)* the shear stress 
distribution found to closely approximate the experimental 
distributions in these inner regions is given by 



This expression is considered to be a closer approximation 
than equation (IV-8) even in the non-asymptotic flow cases, 
and is utilized in the formulation of the sublayer models 
presented in Chapter V. 
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Figure 15. 

U + vs y + , F = 0.002 and K = 0.57 x 10“ 6 , run 42468 
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Figure l6 , 

ft 

U + vs y + , F = 0.002 and K = 1.45 x 10” , run 8l668 
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u + vs y + . 


Figure 17. 


F = 0.004 and K = 0.57 x 10 



run 41268 



Figure l8. 

U + vs y + , F = 0.004 and K = 1.45 x 10' , run 82068 
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Figure 19* 

' f) 

U + vs y + , F = -0.002 and K = 0.57 X 10" , run 52868 
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Figure 20. 

U + vs y + , F = -0.002 and K = 1.45 x 10" 6 , run 80768 
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Figure 21, 


-6 


U/Uoo vs y/6 , F =0.0 and K = 0.57 x 10 3 run 51^68 






Figure 24. 

U/ Uoo vs y/6 , F = 0.004 and K = 0.57 x 10 ^ , run 41268 
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Figure 25. 

0.004 and K= 1.45 * 10 


6 


run 82068 



Figure 26. 

U/Uoo vs y/6 , F « -0.002 and K = O.57 x 10" 6 , run 52868 
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Figure 27. 

U/U* vs y/6 , F = -0.002 and K = 1.45 x 10~ 6 , run 80768 
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Figure 28 . 

Asymptotic values of Reo 0 vs K 







Figure 20. 

Comparison of Simpson's shape factor correlation 
with data 
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0 . 050 X 
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0.0 

-0.00571 

A 
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-0.00473 
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52868-4 

-0.0368 

-O.OO 367 


Figure 31. 

Dimensionless shear stress profiles for K = 0.57 x 10 
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0.1046 

-0.02524 

O 

81668-3 

0.0450 

-0,01561 

* 
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Dimensionless shear stress profiles for K = 1.45 x 10 
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CHAPTER V 


THEORETICAL PREDICTION OF EXPERIMENTAL RESULTS 

The experimental results presented in Chapter IV are 
restricted to near-asymptotic flows over a relatively wide 
range of pressure gradients and blowing fractions. It is 
desirable to correlate these data in such a manner that, 
for a first approximation, extensions can be made to non- 
asymptotic flows including variable boundary conditions at 
the surface and variable fluid properties. A finite-difference 
prediction program, of the Patankar-Spalding type [ 38 ], was 
developed for this purpose in parallel with the present 
experimental study. Two different semi-empirical inputs 
were formulated, resulting in essentially two different 
programs. These appear as differences in the sublayer models 
and the associated eddy-viscosity relations for the outer 
regions of the layer. One sublayer model can be described 
as a two-layer model whereas the other is a variation of the 
damped eddy-viscosity distribution proposed by Van Driest 
[ 39 ]. 

A. Patankar-Spalding Procedure 

The procedure provides for the numerical solution of 
simultaneous parabolic differential equations which describe 
turbulent boundary layers and associated transport phenomena. 
Although the transport of energy, chemical species, etc, are 
considered, attention is concentrated on the transport of 
mass and momentum in the present application. 

The boundary conditions at a solid, surface are imposed 
by means of sublayer correlations which are patched to a 
finite-difference solution. This avoids the numerical dif- 
ficulties inherent in regions of large velocity gradients. 
Conservation of mass is forced to be satisfied in this 
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patching procedure, as well as in the numerical solution 
applied to the outer regions of the layer. 

The numerical procedure is advanced in the streamwise 
direction in an amount which depends on the amount of fluid 
entrained. Modifications of this entrainment procedure are 
available in the event additional transport phenomena is 
considered. 

In the present application only two-dimensional boundary 
layer flows are considered, but the procedure itself can be 
extended to other flow systems. Axially-symmetric external 
flows, internal flows, wall jets, and jet flows are examples 
of the flow systems where the procedure has been used for 
purposes of prediction. 

B. Semi-Empirical Inputs to Prediction Program 

In order to have a complete set of equations to solve, 
empirical rate equations must be used; the specific need in 
the present application is the relation between effective 
shear stress and velocity gradient. Appropriate correlatic s 
of the present experimental data are obtained for this purpose. 
As indicated these physical inputs to the program are presented 
in the form of sublayer models and associated eddy-viscosity 
distributions in the outer regions of the layer. The two 
sublayer models are developed, and the associated eddy- 
viscosity relations utilized in the outer regions are described. 
A summary of the semi-empirical relations obtained is given 
in Appendix C. 

B. 1. Sublayer models 

The turbulent contribution t+ to the effective 
shear stress t is calculated on the basis of Prandtl's 
mixing length hypothesis: 

„. -2.dU|dU .. 

T t = P £ I5yl5y ( V - l} 
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The two models considered are in the form of assumed mixing 
length distributions. 

a. Case of constant free-stream velocity 

Simpson [17] has shown that the fully turbulent 
portion of boundary layers with blowing but no acceleration 
can be fitted by the bi-logarithmic law 

7 ? H 1 + t 1 + n O V2 ] - 07W in lfrl 

" (V-2 

The form of this relation is based upon the shear stress 
distribution given by 


x + = 1 + V + U + (V-3) 

W 

This was used in connection with equation (V-2) in arriving 
at mixing length distributions in the present study. In the 
case of suction representative profiles of Simpson, shown 
in Figure 33* were utilized in place of equation (V-2). 

i. Two- layer model 

In this model, the mixing length is assumed to 
have the following form in the inner regions of the layer 

0 < y + < yj i = 0 

yj < y + & - K y 

Here, the mixing length constant k is taken to be 0.44 in 
order to be consistent with the results of Simpson. 
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The critical thickness y+ at which the profiles can be 

C -f- 

considered fully turbulent is correlated, with V w . This 
was accomplished by an iterative scheme which determined the 
positio/i at which equation (V-2) is satisfied simultaneously 
with the viscous sublayer model. In the case of suction, 
the values of y c were obtained, by matching two representative 
profiles of Simpson as shown in Figure 33. Here, a Runge- 
Kutta procedure was utilized, in the solution of the first 
order differential equation encountered. 

The resulting correlation of y + with V + is fitted 
by the following equation, where y+ is based upon the local 
value of shear stress given by equation (V-3). 

5f + = ll.O - 18.0 V + (V-4) 

c w v 

This simple correlation indicates a decrease in the laminar 
sublayer thickness with blowing and an increase with suction 
(in the non-dimensional sense). 

ii. Continuous modified Van Driest model 

The two-layer model does not properly handle 
that region of the profile where both laminar and turbulent 
viscosities are comparable, due to the inherent discontinuity 
in this region. For accelerated flows which are dominated 
by the inner regions of the layer, this has a significant 
effect on the prediction results. An alternative model was, 
therefore, considered in the present study which has a 
continuous distribution of mixing length. As proposed, by 
Rotta [ 40 ] , the mixing length in transpired boundary layer 
flows is assumed to have the distribution given by 

0 < y + i = xy 



78 



Here, A* is correlated with V w + . This damped mixing 
length was first proposed by Van Driest [38] as a good fit 
to impermeable flat plate data using A* = 26 , 

Using the shear stress distribution given by equation 
(V-3), the required correlation of A* with V w + was obtained 
by matching the resulting velocity profiles in the fully 
turbulent region with those predicted by means of the two- 
layer model. The resulting correlation is given by 


A* 


26.0 - 88.0 V + + 110.0 
w 


26.0 - 88.0 v w + + 210.0 



(V-6) 


The better agreement, with the experimental profiles, obtained 
with this continuous relation relative to that found with the 
two-layer model is shown by the data matching comparison given 
in Figure 32 for the suction data of Simpson [17]. 


b. Effect of acceleration 

In both models, the effect of acceleration is 
assumed to be linear in terms of the pressure gradient 
parameter p + for a given value of V M + . Hence, the assumed 
forms of the desired correlations are given by 


and 



A* 




(V-7) 


(V-8) 
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The functions Qy(V w + ) and Q*(V W + ) were determined from 

the present experimental data. 

Eighteen experimental profiles, representing the char- 
acteristics of the boundary layer flows presently considered, 
were used to determine Qy(V w + ) and Q*( v w + ) . The flows 

in which similarity was not attained in the inner regions 
were not considered. 

As indicated in Chapter IV, the shear stress distribution 
in strongly accelerated flows is best approximated by the 
following equation 



Utilizing this relation and the assumed forms of the Prandtl 
mixing length distributions, the experimental profiles were 
matched in the fully turbulent region by an iterative deter- 
mination of the required values of y + and A* . The experi- 
mental profiles and corresponding distributions in accordance 
with the two-layer model and continuous model are presented 
in Figures 34 through 40. The required values of y+ and 
A* are also presented with each profile. 

The continuous model gives the closest approximation 
to the data in the region near the wall (y + < 70), as antici- 
pated, The experimental data are bracketed by the different 
models in this region. The prediction of the experimental 
profiles with the continuous model is also shown to improve 
with blowing. 

The values of Qy(V w + ) and Q*(V W + ) corresponding to 
these data are presented as a function of V + in Figures 

W 

4l and 42, respectively. The continuous distributions fitted 
to these values are also present in the corresponding figure. 
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The scatter of the profile values about these assumed distribu- 
tions appears to be quite large due to the expanded scale 
utilized. The deviations indicated are still within the 
experimental uncertainty of the data. 

Equations (V-7) and (V-8) are utilized in the present 
prediction of these data, with the functions Qy(V w + ) and 
Q*(V W + ) being given by the continuous distributions indicated 
in Figure 4l and 42, respectively. The shear stress distri- 
bution given by equation (V-9) is also assumed to apply in 
the inner regions of the layer not computed by means of the 
finite-difference procedure. 

These correlations reduce to the zero pressure gradient 
cases for p + = 0 . In the event that V w + = 0 and p + = 0 , 
the flat plate "law of the wall" is obtained with both models 
and associated correlations. Hence, these sublayer relations 
are consistent with accepted turbulent boundary layer behavior 
in the more restrictive flow systems. 

B.2. Eddy-viscosity distributions in outer regions 

The assumed mixing length distributions applied to 
the inner regions of the boundary layer obviously cannot be 
extended to the outer regions. It is, therefore, necessary 
to assume a different form of the mixing length distribution 
in the outer region. 

Escudier [40] found that in a variety of boundary layer 
flows the mixing length distribution can be fitted by the 
relations 


y/6 < A/k 

& - xy 

y/6 > A/k 

i = A6 


(V-10) 


Here, the quantity A is taken to be a constant whose sug- 
gested value is 0.09. A similar truncation of the mixing 
length is utilized in the present predictions. 
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The eddy-viscosity data of Simpson [4l] obtained for 
a wide range of blowing conditions indicates that A should 
not be considered a constant for low values of Reg^ . A 
better approximation is given by the relation 


A = 0.25 Reg^ 0 * 125 (V-ll) 

for Reg 2 < 56 OO.O . For larger values of Reg 2 , it is sug- 
gested that A be truncated at a value of O.O 85 . 

In connection with the two-layer model, the mixing length 
distribution given by equations (V-10) is utilized. The 
parameter A is also taken to be a function of Reg 2 as 
indicated. In addition, it was found necessary to apply a 
blowing correction to equation (V-ll) in order to accurately 
match the theoretical profiles with experimental profiles 
in the prediction of highly blown boundary layer flows. The 
resulting relation utilized for A is given by 


A = 0.25 Re 6 2 0,125 [l - 67 . 5F] (V-12) 

where A is again truncated at a value of O.O 85 . 

In the continuous model the mixing length distribution 
assumed is given by 


y/5 < A/k 




y/6 > A/k i = A6 


(V-13) 


Here, A is also assumed to vary with Reg 2 and. F in 
the manner utilized in the two-layer model. 
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Prediction of Constant Free-Stream Velocity Data of 
Simpson 


r* 

S-/ • 


The present program was used to predict the constant free- 
stream velocity data of Simpson [17] in order to check the 
program. The constant blowing fraction runs considered consist 
of two blowing runs, an impermeable wall run, and a suction 
run. 

In order to demonstrate the similarity obtained between 
experimental and theoretical profiles, the theoretical profiles 
predicted for flow over an impermeable flat plate are presented 
in Figure 43. These U + vs y + profiles, corresponding to the 
two-layer, model are shown to be "normal" and acceptable; 
similar profiles are obtained with the continuous model. 

In Figure 44, the friction factors predicted for all 
zero pressure gradient runs considered are compared with the 
experimental data. The values of Cf/2 predicted by the 
two-layer model and continuous model are shown to agree with 
the experimental values at corresponding values of Reg 2 . 

For the highest blowing run, slight deviations are noted but 
in considering the degree of experimental uncertainty in 
these data the deviations are found to be acceptable. 

It is concluded that the prediction program yields 
acceptable results in zero pressure gradient flows. 

D. Prediction of Present Experimental Results 
D. 1. Initial profile input 

The Patankar-Spalding procedure requires an initial 
velocity profile and initial profiles for additional transported 
quantities considered. For prediction of the present data, 

1/7- th power velocity profiles were assumed for this purpose 
in all cases. The corresponding thickness of the boundary 
layer was chosen such that the predicted momentum thickness 
coincided with that experimentally obtained value at a given 
station in the constant free-stream velocity approach region 
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of the duct. This matching of momentum thickness was made 
in the impermeable wall cases and the corresponding initial 
profile was assumed to be invariant with blowing or suction. 
This approach was considered reliable in that computations 
were begun at the upstream edge of the porous region where 
the layer has yet to adjust to the blowing or sucking condi- 
tion imposed. 

D.2. Input of blowing rates and free-stream velocity 
distributions 

It is necessary to stipulate the boundary conditions 
corresponding to the boundary layer flow to be predicted. 

The asymptotic and near-asymptotic boundary layer flows of 
interest were characterized by constant values of K and 
F chosen to fit the associated experimental distributions. 

The free-stream velocity distribution was determined by 
means of analytical expressions representing the streamwise 
distributions of the parameter K . The initial value of 
the free-stream velocity corresponded to that experimentally 
determined for the constant free-stream velocity approach 
region of the duct. Sinusoidal distributions for K were 
used to proceed from one constant K region to another in 
a continuous manner. 

The mass flux at the wall was described by analytical 
expressions so that the required constant blowing fraction 
was maintained with the associated free-stream velocity dis- 
tribution, 

D.3 . Comparison of prediction with data 

The theoretical predictions of H , R e 62 ’ an< * 

Cf/2. for the present data are compared with the experimentally 
determined values in Figures 45 through 54. These quantities 
are presented, as functions of the distance from the upstream 
edge of the porous section. The assumed distributions of the 
pressure gradient parameter K are presented in all cases. 
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and. for the data runs where K = 1.45 x 10" , the resulting 

free-stream velocity distribution is also compared with the 
experimental values at the given traverse stations. 

With the exception of runs 42468, 52868, and 80768, the 
theoretical predictions for both models are found to be good. 
In all the predictions, with those exceptions noted, the 
trends of the data are distinctly represented and the agree- 
ment with the data in terms of H and C^/2 is within the 
calculated experimental uncertainty of these parameters. The 
associated predictions of Re5 2 are also acceptable -or 
most engineering purposes. 

For run 42468 (where K = 0.57 x 10“^ and F = 0.002), 
the predictions are shown in Figure 47 to underestimate the 
value of Cf/2 at the traverse station in the approach 
region of the duct. The experimental value of C^/2 at 
this value of Regg is found to agree with the corresponding 
Cf/2 data of Simpson [17]; therefore, it is suggested that 
this discrepancy is the result of insufficient adjustment of 
the initial profile. For run 52868 (where K = 0.57 x 10"6 
and F = -0.002), a similar discrepancy is shown to exist in 
Figure 50. Here, the theoretical predictions of Cf/2 are 
found to be in close agreement with the corresponding data 
of Simpson; this discrepancy is not considered to be the 
result of inadequate data correlations. 

Run 80768 (where K = 1.45 x 10"^ and F = -0.002) is 
shown in Figure 54. An adequate prediction of this case was 
not excepted with the present program. The experimental 
profiles are known to have different characteristics than 
those required by the semi-empirical models utilized. It 
is presented here only to indicate the limited range of each 
physical model proposed. The two-layer model would not 
function at all in the presence of such strong pressure 
gradients. The continuous model did allow prediction of 
the experimental data but poor agreement was obtained as 
indicated. 
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An obvious comparison to be made is that of the relative 
merits of the two sublayer models proposed. In the prediction 
of Regg , the models yield similar results whereas deviations 
are observed between the two for H and C^/2 . The continuous 
model yields the closest estimate of shape factor H as 
expected due to its good agreement with the experimental 
profiles. The two-layer model is found to yield tne closest 
estimates of Cf/2 in the pressure gradient flows. Friction 
factors are equally well predicted by both models in the zero 
pressure gradient regions. 

In conclusion, it is found that in terms of an overall 
prediction one model does not yield, superior results relative 
to the other. If either H or C^/2 are of major importance, 
the appropriate model should be utilized. 

The development of Re5 2 is overpredicted in the recovery 
region (constant free-stream velocity) for the strongest 
acceleration. Upon inspection of the prediction results, it 
was found that a momentum unbalance exists in the prediction 
for this region of the flow whereas conservation of momentum 
is satisfied in all other regions. 

E. Summary of Results 

Theoretical predictions of the present data were obtained 
with the Patankar- Spalding procedure [38] using semi-empirical 
eddy-viscosity relations. Two eddy-viscosity models were 
considered using Prandtl's mixing length hypothesis. The 
inner regions of the boundary layer in one are represented 
by a two-layer model and in the other by a continuous modified 
Van Driest model. The outer regions are represented, by means 
of a truncation of the mixing length. The effects of transpira- 
tion and acceleration are incorporated in terms of the parameters 
V + and p + , respectively, in the inner regions. The 
bi-logarithmic "law of the wall" proposed by Simpson [17], 
the suction data of Simpson, and. the present asymptotic and 
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near-asymptotic boundary layer flow data were used to formulate 
these models. 

Theoretical predictions of experimental constant blowing 
fraction runs of Simpson in zero pressure gradients are found 
to be acceptable. 

The theoretical predictions of the present pressure 
gradient data are similar ily found to be acceptable, except 
in those cases where the turbulent models were known a priori 
not to apply according to the present data. It follows that 
the proposed models are reliable for prediction purposes 
in these types of flows. 

In terms of an overall prediction both models appear 
equally useful. The continuous model yields the best estimates 
of shape factor whereas the two-layer model yields the best 
estimates of friction factor. 

Some physical implications of the successful prediction 
of these data with such simple models are noteworthy Implicit 
in these models are the effects of transpiration and pressure 
gradients on the increase or decrease of eddy-viscosity near 
the wall. Although these effects are coupled, the same 
qualitative behavior is found in those cases where either 
transpiration or acceleration are considered separately. 

In addition, these are turbulent models which suggests that 
the present flows can be considers i turbulent in nature, 
except in those cases of suction, noted in Chapter IV, where 
substantial structural changes are suggested in the inner 
regions of the layer. 




Figure 33. 

Velocity profiles for cases of suction at wall 



Figure 3^. 

Velocity profiles for cases of suction at wall 
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Figure 35. 

Velocity profiles for impermeable wall cases 



Figure 36 . 

Velocity profiles for impermeable wall cases 
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Figure 37 . 

/elocity profiles for cases of blowing at wall 



Figure 38 . 

Velocity profiles for cases of blowing at wall 
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Figure 39- 

Velocity profiles for cases of blowing at wall 



Figure 40. 

Velocity profiles for cases of blowing at wall 
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Figure 4l. 

Correlation of Q with V + 

jr W 



Figure 42. 

Correlation of Q* with V w + 
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Figure 43. 


Predicted U + vs y" 1 profiles for impermeable 
flat plate utilizing 2-Layer model 
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Figure 44. 

Prediction of C f /2 vs Regp according to data 
of Simpson 
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Figure 51. Comparison of prediction with data: run 8l668, F = 0.002 
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CHAPTER VI 
SUMMARY 

One-hundred- twenty-six boundary layer mean velocity 
profiles from the Stanford Heat and Mass Transfer Apparatus 
are reported for accelerated flows with injection or suction 
at the wall. The pressure gradients investigated are those 
where K = 0.57 x 10" 6 , 0.77 x 10" 6 , and 1.45 x 10~ 6 . 

For each pressure gradient, the boundary conditions investi- 
gated cover a range of constant blowing and sucking fractions 
from F = -0.004 to 0.006 . 

These data serve as the basis for eddy-viscosity models 
which, combined with a finite-difference procedure, result 
in predictions of the data. 

A. Conclusions 

A . 1 . Experimental data 

a) "Best estimates" of friction factors, based on 
the mean velocity profile data, are found to be self-consistent 
and in acceptable agreement with other experimental deter- 
minations, according to either the momentum integral method 
or the viscous sublayer model method. It is found that the 
present constant free-stream velocity friction factor data 
agrees with that of Simpson [17]. The effects of acceleration 
are found to be different in the cases of blowing and suction, 
but involves only small deviations from the following cor- 
relation of Simpson. 

Q yf = 0.0130 [ ^ n -ll +5 .l-] °‘ 7 (R e 5 2 ) ^ 

In suction, the friction factors decrease, whereas with 
blowing an increase is indicated. 
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b) in the impremeabie wall cases, similar profiles 

in the inner regions of the layer were attained for all three 
values of the pressure gradien' .orameter K which were 
investigated. Two characteris & of the boundary layer are 
shown: (l) The profiles deviate from the flat plate "law of 

the wall" in an "overshot" manner within the logarithmic 
region, and (2) the wake region is substantially diminished. 

The degree of the profile "overshoot" in the logarithmic 
region and the associated decrease in the wake was found to 
be greater for progressively higher values of K . The 
extent of the logarithmic region, in U + vs y + coordinates, 
is also found to decrease as K increases. 

These characteristics may be the result of surpressed 
turbulence in the region near the wall and a reduction of 
shear stress in the wake region, although turbulence measure- 
ments were not taken. 

c) The qualitative characteristics of the inner 
regions for boundary layer flows in the presence of a strongly 
accelerated mainstream apply for the range of blowing and 
sucking fractions considered. The quantitative effects are, 
however, dependent upon blowing or sucking fraction. Similar- 
profiles are found to exist in the inner regions, except in 
flows where the structure of the layer appears to be sub- 
stantially changed, i.e., where the layer may be re laminar izing. 
In the majority of the suction runs, non-similar behavior of 
this latter nature was observed. 

d) In the outer regions of the boundary layer, 
similar profiles were attained for all bloving and sucking 
fractions considered, with the exception of the flows where 
the sucking fraction was F = -0.004, In this exceptional 
case, the profiles are found to continually adjust toward a 
seemingly laminar mode. It is noteworthy that similarity 
is also implied in "velocity-defect coordinates" as well as 
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U/Uqo vs y/6 . This is a result of the friction factor being 
essentially constant in the flow direction. 

This behavior in the outer regions coupled with a similar 
conclusion relating tc the inner regions confirms the existence 
of completely similar profiles in asymptotic boundary layer 
flows with the exception of those cases noted. 

e) The boundary layers reported here are believed 
to be turbulent. Unique values of Regg and shape factor H 
correspond to each set of K and F for asymptotic boundary 
layers. For asymptotic flows on an impermeable wall, both 
Re &2 and H are found to be in better agreement with those 
values predicted on the basis of an assumed turbulent 1/7- th 
power profile than with those for an exact laminar solution. 

f) The shape factor correlation of Simpson over- 
predicts the experimental data and is, therefore, not considered 
applicable in such strong favorable pressure gradient flows. 

g) Shear stress distributions were calculated for 
the present mean velocity profile data in the asymptotic and 
near-asymptotic flow regimes. It i’ found that the maximum 
shear stress occurs for values of y+ less than 25 . 

In asymptotic or near-asymptotic boundary layers, the 
shear stress distribution in the vicinity of the wall is found 
to be approximated by the dimensionless relation 



In the event p + = 0 , this reduces to the shear stress dis- 
tribution assumed in the formulation of Simpson's bi-logarithmic 
"law of the wall". 
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A. 2. Theoretical prediction 

a) Either of two simple Prandtl mixing length models 
formulated predict profiles in acceptable agreement with those 
experimentally obtained. 

b) The proposed models, summarized in Appendix C, 
were used in a Patankar-Spalding prediction procedure. 

Predictions of H , Re& 2 , and C f /2 for the present data 
are found acceptable for most engineering purposes. 

c) In term of an overall prediction, neither model 
appears superior to the other. The continuous modified Van 
Priest model is observed to yield the best estimates of shape 
factor H whereas the two-layer model yields the best estimates 
of friction factor. 

B. Recommendations for Future Work 

B. l. Experimental research 

a) Execute a series of asymptotic runs at suf- 

fi iently high accelerations to include complete relaminarization. 

b) Non-asymptotic flews should be investigated which 
include variable accelerations and wall boundary conditions to 
determine the restrictions on the results of the present study, 

c) The decelerated boundary layer with blowing and 
suction should be investigated in order to complete the present 
study which was restricted to accelerated boundary layers. 

d) Hot-wire anemometer studies of the fluctuation 
components of velocity should be undertaken, to study the 
structure of the layer resulting from these coupled effects 
of acceleration and transpiration. This investigation would 
also allow the generation of shear stress profiles with which 

to compare those generated in this work and that of Simpson [IT]. 

e) Foreign gas injection poses no problem for the 
operation of the present apparatus. Suitable techniques are 
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available for sampling and measuring concentration profiles. 
Such data would materially aid in studying the relati/e 
effects of molecular and turbulent transport mechanisms. 

B.2. Theoretical research 

a) The present prediction model should be modified 
to handle highly non-asymptotic characteristics on the basis 
of data obtained in the proposed research investigations. 

b) Utilizing the present prediction program, a 
systematic investigation of the influence of high velocity 
and large temperature difference should be undertaken. 

c) Other existing prediction procedures should be 
combined with the required correlations of the piesent data, 
if possible, and resulting predictions of the present data 
should be obtained. Upon comparison of these results with 
the present predictions, the most appropriate model could 

be chosen. 


Author's Note: The programs mentioned in items 

B.l. a,b,d and B.2. a,b are presently 
being pursued. 



APPENDIX A 


TABULATION OF EXPERIMENTAL DATA 

The experimental data is presented in the following 
order according to blowing fraction: 


I. Unblown flat plate runs 


U<« (ft/sec ) 
42 
42 
86 
126 


Date of run 
31067 
71967 
72067 
80867 


II. Pressure gradient runs 


Date of run 


F 

t 

K = 0.57 x 10" 6 

0.77 x 10 -6 

A 

1.^5 X 10" 6 

0.0 

51468 

51568 

73068 

0.001 

51068 

50768 

81368 

C .002 

42468 

42368 

81568 

0.004 

41268 

40268 

82068 

0.006 

— 

120867 

82668,82268 

-0.001 

52068 

52168 

80568 

-0.002 

52868 

52768 

80768 

-0.004 

60168 

60468 

80968 
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For each blowing fraction , the experimental data is 
presented in the following order: 

SUMMARY TABLE 

This table contains the friction factors, integral 
parameters, blowing or sucking fractions, and conditions 
relevant to each profile. Uncertainty estimates for these 
quantities are presented for each run. 

The enthalpy thickness discrepancies are also presented 
where the subscripts p and St refer to profile and 
Stanton number estimates, respectively. 

SETUP DATA 

Ambient conditions, free-stream velocity distributions, 
and A" distributions for each pressure gradient run are 
tabulated. 


VELOCITY AND SHEAR PROFILES 

Velocity and calculated shear stress profiles are 
tabulated. 
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SWAT FOR 

THE FOLLOWING PROGRAM REDl'CFS CONSTANT VELOCITY PROFILE OATA TAKEN 
ON THE STANFORD HEAT AND WAS TRANSFER APPARATUS 


SOME OF THE SYMBOLS UTILIZED IN THE PROGRAM REPRESENT THE FOLLOWING 
IMPORTANT QUANTITIES l 

CFFLAG- l.C IF PROFILE IS IN PRFSSUAE GRADIENT REGION 
0.0 IF NOT IN PRES5URF GRADIENT REGION 
CHFLAG* 1.0 IF LARGE ROTAMFTSPS U5E0 IN FORWARC SECTION 
0.0 IF SHALL ROTAMETERS USED IN FORWARD SECTION 
NPLATE* PLATE AT WHICH DIFFERENT ROTAMETER IS FIRST USED 
0 IF ONLY ONE TYPE OF ROTAMETER USED 

CM- ROTAMETER 5ETTINGS 

EO- PLATE TEMPERATURE (MILLIVOLTS) 

00- HEIGHT OF PROBE HEAD ( INCHES ) 

P- TOTAL PRESSURE MINUS REFERENCE STATIC PRESSURE (IN. H20I 
NOTE: REFERENCE STATIC PORT ADJACENT TO PROBE TIP 
PBAR- BAROMETRIC PRESSURE (IN. HGI 
PQ- REFERENCE STATIC PRESSURE (IN. H20I 

PRGT- PRESSURE OF ROTAMETERS ABOVE AMBIENT PRESSURE (MM HG) 

PTOTAL- TOTAL PRESSURE IN DUCT (] N. H20) 

RHUM- RELATIVE HUMIDITY 

TAMS- AMBIENT TEMPERATURE (OEG-PI 

TG* FREESTREAM GAS TEHPEAATURF (MILLIVOLTS OR DEG-F1 
TROT- TEMPERATURE OF GAS ENTERING RDTAMETTP 5 (MILLIVOLTS) 

X- DISTANCE DF PROBE TIF FROM ENTRANCE ''S 1 E5T SECTION (INCHES) 

V* DISTANCE FROM PLATE SURFACE (INCHES) 

U- MEAN VELOCITY (FT/SEC I 

UGG- LOCAL FREESTREAM VFLOCITY (FT/SEC) 

UUG- U/UGG 
YOEL- Y/OEL 

OEL* Y-P0S1 TION WHERE UUG-C.M (INCHES) 

0ELT1- DISPLACEMENT THICKNESS (INCHES) 

0ELT2- MOMENTUM THICKNESS (INCHES) 

H> SHAPE FACTOR 

REQEL2* MOMENTUM THICKNESS REYNOLDS NUMBER 
PM- LOCAL BLOWING FRACTION VO/UGG 

MOOTH- LOCAL HASS FLUX AT PLATE SURFACE ( LBM/SEC-FT2) 

KM- LOCAL PRESSURE GRADIENT PARAMETER 

CF2SL- CF/2 BASED ON SUBLAYER EQUATION AND FIRST OATA POINT 

CF2SL2- CF/2 BASED ON SUBLAYER EQUATION AND SECOND OATA POINT 

UPLUS1- U* BASED ON CF2SL 

YPLUS1- Y. BASED ON CF2SI 

PPLUS1- P* BASED OH CF2SL 

VPLU51- VO* BASED ON CF2SL 

CFMICR- CF/2 BASED ON COMPLETE TRANSFORMED MOMENTUM INTEGRAL EQU. 

CF2MI- CF/2 BASED ON ASYMPTOTIC FORM OF MOMENTUM INTEGRAL EQUATION 

UPLUS2- U* BASED ON CFMICR 

YPLUS2- U* BASFii ON CFMICR 

PPLUS2- P* BASED ON CFMICR 

VPLUS2- VO* BASED CN CFMICR 

1. DECLARATION OF VARIABLE TYPES 

REAL At BtCiCFCORRI 10 )i CFFLAGflO ( t CFMICRI10) (CF2MI ( 10) ,C F2SLI 10 1 . 
ICF2SL2I 10) tCH20,CH(2A) ,CMFLG,CPl,CP2,D,0CP,DDP 1 10) , DEE1 (10 1 , 

2DEE21 10 1.DEl ( 10 ) iDEL 1( 10 1 ,0EL2( 101 ,0ELT1( 10) .0ELT21 101 , DUDXIAB ) , 
3E0(2AI . EPS , F 1 24) .FM( 10 I , H(10) *K(4B) ,KFL0W(24) iKFUDGEI 24) *KM{ 10) t 
AMAtMDOT I 24 1 , MOO TNI 10 I ,MV,OD< 10) ,P(10.501 . PAMB tPBAR.FI NF (10) tPO(AB) 
5 tPPLUSlt 10 1 .PPLUS2 (10 ),PROT,PPOTA. PROTAB, PTOTAL »PSAT( 12) ,PVAP,RA , 
6RAD,REDEL2tl0l,REP,REPl,REPS.REX(A8),REXUGGtl0 ),RH0A,RH0GIA8 (,RHOH 
7 ,RHOL t RHOSAT 1 12 ) t RHOUt RHOVtRHOIROI 24) »RHUM»RN, SFLAGI 10) , 
BSUMlIlOtSOt <SUH2(10t501 ,T ,TAMB, TAVGI2A ), TEMPI 12 ) ,TG,TO( 2A I .TROT, 
9TR0TA.UI 10.50) ,UGf 4SI ,UGGI 10 ) .UPLUSI 1 10 .50 1 .UPLUS2 (10 ,50 ) , 

1 UUG ( 10,50 ) . VAFH.VAPL.VEP5. VISCI 10) ,VISCG( AB ). VPLUS1 (101 , VPLUS2(10t 
2,VZER0(2A) , YZERCNt 10 1 ,WAC T{ 24 > ,WSTD(2A ) . WSTDJ. X 1 10 1 , XS ( AB) « 
3Y(10,50) tYDELI 10,50 ) , Y0( 10) .YPLUS1I 10. 50) . YPLUS2 ( 10,50) 

INTEGER OATE.DP.NPREF (10 ) .R.RUNNO.S. 2110) 
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II* FIXED OAT 4 INPUTS 

A: KFLUM , KFUOGE* CORRECTION TEAMS THAT WEIGHT ROTA RE T£ R READINGS F OR 
FLOW THROUGH CFNTFR 6- IN. SPAN OF PLATE 

DATA KFL0WI1I .KFLOWI 2),KFL0W!3I ,KFLQW(41 , KFLOWI 5 1 i K F L OH ( 4 I , KFL OH ( 7 
II ,KFL0wm,KFL0WI9) .KFLOWtlft ),KFtOW( 11 1 ,KFL0WI 12) ,K E LQWmi , 
2KFL0WI14I .KFLOWU5) ,KFLOWIl6l .KFLONUTI .KFLOWI 1 91 ,KFL0WI19), 
3KFLOH) 20 I .KFLOWI 21 1 ,KFLOHI 22 1 »KF LOW (211 .KFLOWI 241/1* 020 A ,l«Qiri, 
41.0309, 1.3417,1. 0309, l.C 309, 1.0183,1. 0493,1 .C3»5, l. P449, 1.0331, 
51.042 8, 1.05C4, 1.0373, l.f 526, 1.0152, 1.0341, 1 .C331 ,1.0081 ,i.t 471 , 
61.0363, 1.0428, 1.0018,1.0331/ 

DATA KFJDGc II I .KFUOGE 12 1 • KFUOGE ( 3 I .KFUDGE (41 .KFUDGE I 51 tKFUOGEI 61 , 
1KFU0GE I 7I.KFUDGEI 8 > .KFL'DGEI 91 .KFUDGE 1 1 0 > , KFUOGE ( 1 1 I ,KFUDGEI 121 , 
2KFU0GEI 1 3 1 , KFUOGEI 14 1,F FUDGE ( 15 ) ,KFU9GE I 161 ,K FUDGE I 1 7 I ,KFUOGE IIP!, 
3KFU0GEI 19* .KFUDGE ( 20 1 «► FUDGE! 21 1 , KFUOGE (22 t , KFUOGE I ?31, KFUDGEI 2‘ > / 
4 -*I),U IP ,0.924*0.11,-0. CC/ 5 , 0.008C , ft ,304, 0. 004 , -0. 90S , ft. 008,0.0,0 .I" 99 

5. 0. U08, )..), u. UU, 0.006,0.016, 0.010,0.01 6, 0.016,0. OP 3, P, 016,0.01”, 
60.010 ,U*U08 / 

AS XS' POSITIONS OF STATIC PCRTS 

DATA XSIll,K$m,XSI3>,XSI4),XS(5l,KSt6),X$(TI,X$tRI,XSf9),XStlPt» 
1XSI11 l,XSI12l ,XS( 13I,XS(14I,XS(1S* ,XSI16* ,XSI 1 7 I ,XSU8 I ,X5I19I, 
2XSt20>.X$I21),XSt22l»XS!23>.XSI 24) , XSI 2 5) , XSI 26 1 , X5I 271 ,KSI 281 , 
3XSIZ9l,XSt3C)»XS(31l»XSI32l,XSI33l»XSI34l,XS(35l,XSI36l,XSI37l, 
4XS( 38 1 , XS( 391 • XSI40* ,XS! 41 1 , XSI 42 1 , XSt 43 1 , XS I 44 * , XJ 1 45 ) , XS 1 46 1 , 
5XSI47 1 , XSI48 )/l .969, 3.933 3.933,7.961,9.969,11.953,13.937,13.943, 
617. 933,19. 922,21.938,23.954,25.962,27.962,24.978,31.939,33.953, 
735. 95S, 37.971,39. 987, 41. 963, 43. 963, 45. 963, 47. 9 79, 49. 979, 5 1.970, 

853.995,55.971,57.971,59.955,61.479,63,971,65.979,67.963,69.971, 
971.979,73.963,75.934,77.947,79.439,51, 931 ,83. 962 ,A5. 931 ,07.915, 
289. 939, 91. 931, 93.947, 96.CZ 

Cl PROPERTIES OF ORY UR 

DATA TE3imi.TEHPI2l.TEMF (3 1 ,TEHP* 4) ;TEHP( 5 t ,TEHP (6 ( , TEMP IT) , 
ITEMPtSI .TEMPI 9 1 /40.0,5a. 3,60.0, 7ft. 0,80. 0,90. 0,100. ”,110.0,1 2C.0Z 
DATA PSATI1 )»P$AT 1 2 l.PSAT 13 1 »P5AT I 41 ,PSAT 1 5) ,PSAT( 6) ,PSAT1 7 * , 
1PSATI8I ,PSAT| 91/1 7.53,25. 65,36. 90, 52. 2ft, 73.00,100.40.136.50,183, 60 
2,243.70/ 

DATA RHOSATIl I , RHOSATI 2 > , RHOSAT 131 .RHOSAT 14 I ,RHOSAT I 5 ), R HO SAT I 6* , 
1RH0SATI Tl. RH0SATI8*,RH0SATI 91/0.000409, 0.000587, 0.000 B30, 0.001153, 

20.001 58 3. 0. 002 13 9.0.00285 3,(1 .00 3 TTQ, 3.004920/ 

01 THERMOCOUPLE CALIBRATION CONSTANTS 

DATA A. 8, C.D/-2220 * 703,781. 2S .7. 950782,0.256/ 

III. READING AND RECORDING OF INPUT DATA 
A: RUN CLASSIFICATION 

1 FORMAT! 16, 2X, II, 2X, II, 2X.F5.0, 63,121 

2 READ! 5, II OATE, RUNNO, S .RAD.NPLATE 

3 FORMAT I 1H1 I 
WRITE! 6,31 
WRITE (6, 3 I 
WRITE I 6, 31 

4 FORMAT 1 3TH RUN NO. NO. OF TRAY. RADIUS) 

HRITEI4.4I 

5 FORMAT I5X,I6,1H~,I1«0X,I1 , IOX «F5, 3) 

WRITE 16,5 ) OATE, RUNNO, S, RAO 
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B ■ I SOI V I DUAL PROFILE ORTA 
DO 13 M*l i S 

6 FORMAT! F6.0, 2K,F6.0.2X, 12, 2X,!2,2X,F5.0, 2X.F4.il 
REA0(5,6I XIMI,PINFIHI,NPAEF(N),2(H),00(MI,CFFLAG(M) 

7 FORMAT I ///,44H X PINF REF l DO CFFLAG1 

WRITE (6,71 

4 FaRHAT<5X,F6.3,2X,F6.4,2X,IZ,3X,l2,2X,f5.3,4X,F4.1,//l 
WRI TE I 6 . 8 I XMl.PINFINI , NPREF ( M I , 1 1 Ml . ODI H I , CFFLAGI M I 

10 FORMAT I 5 X, 21H ft V PT-POREF,/) 

WRITE! 6.10 t 

OP-KMI 
DO 12 R-1.0P 
9 FORMATI 2F1C.0I 

READ! 3.91 Y ( M, R I , P ( M , R | 

11 FORMAT I6X,I2,3X«F5.3,3X,F6.4) 

12 WRITE (6.111 R.YCI.R) .PIM.RI 

13 CONTINUE 
WRITE (6.31 

Ct STATIC PRESSURE DISTRIBUTION DATA 

14 FORMAT 13X.32H I PO(I) I P0( II I 

WRITE 16.141 

DO 17 J-1,24 
L*J*24 

19 FORMAT) 5X. FB.P.24X , F7.G ) 

READ! 5. LSI POIJI.POIL) 

16 FORMAT (5X.I2.3X.F7.4iSX.12.3X.F7.4l 
WRITE (6.161 J.POI J).L.PO(L) 

17 CONTINUE 

Dt SETUP DATA 

18 FORMAT ( 7F8. 0 . F4. 0 ) 

READ! 5. 181 TAM8, TG, TROT, PBIR.RMUM.PROT.PTOTAL.CMFlG 

19 FORMAT if / .5X.55HTAMB TG TROT PBAR RHUM PROT PTOTAl 

1CHFLG) 

WRITE 16 .19) 

20 F0RMAT(5X,F5.2.F7.3.2X,F5.2.2X,F5.2t4X.F4.2,lX,F9.r,lX,F7.4,2X,F4. 
Ill 

WRITE (6,20) TAMS, TO, TROT, P8AR, RHUM, PROT, PTOTAL, CHFLG 

21 FORMATI//, 5X.19H I ECU) CMU)I 
MR ITE (6. 21 I 

DO 24 t-1,24 

22 FORMATI 3X.2F10.0) 

READ! 5,22 ) EOdl.CNIIl 

23 FORMAT (5X f I2,4X,F5.3,3X.F5.2) 

WRITE (6 ,23 ) I.EQdl.CMU) 

24 CONTINUE 

IV. CALCULATION DF FLUID PROPERTIES AND ASSOCIATED CONSTANTS 

At CONVERSION OF THERMOCOUPLE DATA 

1 F( TG.LT, 10.0 I TG-A*e*SQRT(C«C*TG)*49.93 
IF(TROT,LT. 10.01 TROT»A+fl*SORT(C+D*TROT )«49.93 
TROTA«TROT *460.0 

Bl CH20 CORRECTS DENSITY OF H20 FROM 62.4266 FOR AMBIENT TEMPERATURE 

CH20-0, 99732-0.0001395*1 TAM8-75.0) 

Ci PRESSURE CONVERSIONS 

PAHB-PBARP2116.0/29.96 

PROTA-PBAR*PROT/25.4 

PR0TAB»2116.0*PROTA/29.96 
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£ 

C Dt MIXTURE COMPOSITION IS DETERMINED FROM RELATIVE HUMIDITY AND UFEO 
C TO GET MIXTURE GAS CONSTANT RM VIA PERFECT GAS ASSUMPTION 

c- 

DO 25 N-i,9 

I F< TEMP (N| (GEtTAMG I GO TO 26 

25 CONTINUE 

26 T- TEMPI Nl 
EPS-T-TAM6 
VAPH" PS AT ( N I 
VAPL-PSATIN-U 
VEP$*VA PH- V A PL 
RHOH-RHOSATI N1 
RHOL-RHQSATIN-l) 

REPS-RHOH-ftMCL 

RMOVftHOLm0.0-EP$>*REP5/L0.0 

RA-53.3 

PVAP»RHUM*(VAPL*iO.O-EPS J*VEPS/10.0 

RHDA-I ( PAH0-PVAP I / (R A* ( TA PB + 460«0 ) I ♦< RHUH4RH0V ) I 

MV*RHUN*ltHGV/RHOA 

MA-UO-MV 

RM»1545.0*<MA/28,9tMV/ie.GI 

c 

C V, PROCESSING CF FftEESTREAM CATA 
C 

DO 27 I "1 * AS 

RhOGI I A * I °AHS*5 # 20*P0t I H / | RHMTG*46Q,OJ ) 

VISCGf n-(lI.O*O.CIT5*TG ) /< ICOOOOO.GARHOCU 1 I 

27 UG( n*$QRT( <64 # 34*(PT0TAL-P0( m*<62»4*CH20/RK0GU I I/12.0H 
OUDK 4 1 ) ■ < UG ( 2 1 ~UGl 111/1x5(11/12*01 
K(ii»viscGatMouoxm/iuGUi*uG(in 

RE XU 1 *UG< ll*xsa* m2.C*VI5CG<in 
DO 2B I *? * 4 7 

DUOX4 I . =12.0*(UG( 1*1 IHIGI l-LM/tXSI 1*U~X5I MM 

mi u*viscgi n*ouoxin/(UG( vi*ugu it 

28 REXtI1-JG( I )*XS(n/tl2.0*VlSCG (111 

c 

C VI* DETERMINATION OF PLATE TEMPERATURES AND MASS INJECTION RATES 

c 

DO 34 J*l*24 

TGU)»A*e*SGRTIC*O*E0( Jll 

TAVGI JJ *TOI J I *49,97-12, 6E-04AT0< J I- 32, 0E-O6*TO (J » *T0| Jl 
MOOT I Jl-0,0 

IFlCMtJ I .1 E #0 .0 > GO TO 04 
IF ( J. EQ.NPLATE I CMFLG-1 .O-CMFLG 
IF(CMFLG«LE»0«QI GO TO 29 
C 

C LARGE ROTAMETERS J NEW FIT FCR FACTORY CALIBRATION PLUS/MINUS 0.3V 
WSTO J* I 0*60+0* 752*CM ( J I -0 * 50*S I NIC M( Jl* 3* 14 17/25, 01 )*0, 075 / 60. 0 
C 

GO TO 30 
C 

C SMALL ROTAMETERS i NEW FIT FOR FACTORY CALIBRATION PUJS/HJNUS 0,3V 

29 WSTOJ*(0,175*0*13091*CM( J I -0, 06?*$ IN U CN1 J i-2 ,0 1*3. 1417/21 ,Q H* 
10*075/60*0 

C 

30 WSTDI J1-WSTDJ 
C 

C ROTAMETER FLOW CORRECTED FOR DENSITY TO YIELD ACTUAL FLOW* THEN 
C CORRECTED FOR PLATE POROSITY VARIATION 
C 

WACT< JI-WST0UI*SORTHPROTA8/(RM*TRDTA*O,O7SlI| 

I F ( PROT .GT,— 0,01 I MOOT t J 1 *WAC T t Jl * t XFLOW I J 1 * K FUDGE I J I )42, 01230 
I F ( P ROT , LE, -0 * 0 1 1 MOOTUI»UACT«JI*fKFL0WUH#2.Oi23i 
C 

34 JI-2*J-1 

RHQZROI Jl PAMfl+5* 20+P0IJ I 1 I / (AH* ( TAVG( J 1 4460,01 I 
FIJI -MOOT I Jl/IUGI JE »*RHOC< Jill 
34 V I ERO I J I "MOOT ( J I / RHO ZRO ( J I 



c 

C VII. CALC UL A TI CN OP PROFILE VELOCITIES AND ASSOCIATED PARAMETERS 
C 

DO 3B M-lpS 
NP-NPRfcf IM P 

If IKC MI.GE.XSfNPH OOP (Hp — (PdfNP+l l-PO^NP II *( X|«l- X$(NP 1 1 / 
XU5IHP*U-XS(MPM 

IF<X(M)«LT.X£<NPH OOP 1 M I ■ I PQ< NP I-PO( NP-1 U *1 1 $ I NP I - XI H PP / 

1 i XSINPI-XSl NP-1 ) I 
YDt M J »Y I H* 1 P 

IFUfMI.CE* )C< I rjPl) RHOU-RHOGINP !♦< RHQGt NP + 1 1-PHGGINP P >• I MMh-XSl UP 
1 1 1 /I X 5 ( NP+ 1 1 -XS ( NP I I 

I F ( XlHI.LT, KSINPH RHOU«RHQG(NP-lP+|RHOGt NP I -RHQG (NP-L 11*1 K ( M I - 
1KSCNP-1 P )/<XSfNP)-X$INP-l 1 1 

UGG(MP*SQRT(C64, 34*(P1NF (KWOTPCMI I* 1 62 *4*CH20/RHGU I /12 .0 I \ 

I F I X ( M ) . GE » X5 ( NP ) I KK(Mj-MNPl + fMNP+i P-K ( NP P f * f Xf M I -*S < NP U / 
mS(NP*l>-XS(NPII 

IFMtHJ,LT»X$(NPI J KHCHI »*f NP-l I * I * I **P I -K f NP-1 P 1 * I X ! H I -XS ( NP-1 IP/ 

1 t MS tlslp I -KS 4 NP-L I I 
DO 600 J I * 1 1 2 4 
XJI*J I 

XPLD-4,0fXJt 

IF(X(M|,LT.XPLDI GO TO 6G1 

600 CONTINUE 

601 MDOTH< M } "MOOT ( J f I 

V Z E ROM Pi)* V ZERO I J H 
FN(M|*M3CTNIMI/tRH0U*UGGim 

c 

C At ABSOLUTE VELOCITIES AND ASSOCIATED PARAMETERS 
C 

OP-ZIMI 

DO 33 R-l.DP 

P(P,RI»P|N t RI+DOP(MI 

U(tf t R|*$QRTH64.34*P(M*R?*C62.4*CH2Q/ftHQUi/12«rH 

C 

C TOTAL HEAD READINGS CORRECTED FOR VISCOUS EFFECTS 
C 

VlSC(MH(H.040<0175*TGI/UnoOOOO*0*RHOUI 
RFP*UlM,R|*fvAD/l 12*D*V1SC IMP I 
CP2-1.000 

IFIREP.GT. 100*01 GC TQ 61 
REPl-AEP 

IFIPfM f Rl.ie*O.OCUl GO TO 6 L 

31 C Pl*l • 020-0. 566 /SORT C RE PI I *3. 530/ftEPl 
REP1*REP/SQRT(CP1I 

CP2»1 *020-0 * 566/S GRT < RE PI P *3# 530 /ft IP 1 

REP l*ft|P /SORT (CP21 

DCP-CPZ-CP1 

IFIDCP.G1. 0.011 GO TQ 31 
61 U<MfR>»UlMtR)/S0RT(CP2t 
C 

UUG f M * R I *U ( M t R I /UGG f M I 
YlMtftl*YlMtRI-Y0lM|*ODIH|/2.O 
C 

C DISPLACEMENT i MOMENTUM, AND BOUNDARY LAYER THICKNESSES CALCULATED 
C 

IFIR.LE.1P GO TO 40 

|F{UUG<M l R]'GE,0.97*AN0.UUG(NiftJ.LE*O.?9J DEL 1 1 M 1-Y <M , R I 
If (UUGtMiR I #GE •0*97*AND*UUGI H t ft 1 »LE.Q* 991 DE E 1 ( M P ■ 1 * 0- UUG ( M * R I 
1 F t UUG 4 MtR I *GE« 0* 99*AND*UUG (M |R-1 1*LE* 0*991 DEL 2(N I «Y(M»R I 
If<UUG<N ( A*,Gf ,0.99*AN0.UUGfH,ft-n.tE.O.99J 0EE2 (Ml *l*D-UUG iHflH 
1FIR.GT.11 GO TO 32 
40 SUH1 1 MtR P *UUG ( M,R16Y(M|R|/i*0 

$UH2(MifU-UUG<H,RI*UUG(H f R|4YCH t R1/2.0 
GO TO 33 

32 SUM 1 ( M i R I ■ SUM 1 ( M t R- 1 1 ♦ ( UlPG ( M * 1 1 +UU G ( N * R-l H * l Y ( N t R I -V l*s * R-l 1 1 / 2* 0 
SUM2 1 M f R I *SUM2 ( M , R-l 1 + ( UUG( M * ft I HJUG( M < R P + UUG( M, 1 *UUG I H ,R-11 P* 

i(VIMjRI-Y(M ,R— l 11/2.0 

33 CONTINUE 

DEL(M|-DELL<MU(DEL2(MP-OEUlMM*(OEEl(MP-O.DIOP/fDEEliMP-DEE2fM)p 
DELT1 tM|>Y( H t DP l-SUMK M* DP} 

DELT2 < M I -SUM! i M , OP I -SUM 21 M f DP I 
MIM|»Of LT1 IMP/DELT2 (Ml 
C 

REDEL2(MI-UGGlMltOELT2iMI/(VlSClM 1*12.0 I 
REXUGGI M P«UGU(MP*X1M> /( VI SC(M 1*12.0 P 
DO 90 R-t*OP 

90 TOIL (Mi R1«Y IMf RP/OEUMI 
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Bi DETERMINATION OP CF/2 
CALCULATION OF CFZ M 1 
CF2MKMIO.C 

IFICFFLAGIMI.GT.O.OI CF;;M||N|.R£0EL2INI*(HIH1*1.0I*KHIM>-FMM1 

CALCULATION OF CF2SL ANC CF2SL2 

IFINODTMIMI .LT.O. 0001. AND, “DQTMINI.GT. -0.0001 1 GO TO It 
CF2SUN>»KN(NI/FN(M14(FH(M)*UUGIMtU-RN(NI*REXUCGINI*Y|N,ll/X(NI I 
l/(EXPlFHINI*AEXUGGIMI*VIF,lI/XINn-l.OI 
CF2SL2I MI*KN(HI/FMIH lt( FMHI»UUG(M >2 l-KHI H I *RE XUGG I H I *Y IN . 2 I /XI “I I 
l/l£XPIFi(HI*REXUGG(N)*Y<Mt2l/X(M)l-U01 
GO TO 36 

35 CF2$UM»«UUG(M,1 l/(REXUGG(N)*Y(M,l l/XIM) l*0.5*KN(NI*REXUGG(NI* 
1YIK,1)/X<NI 

CF2SL2INI-UUG(M.2I/IREXUGG(NI*Y(M,2)/XINI I tO, 5**HNI ME XUGGIH!* 

1YI H.2 1 /X1NI 

36 CONTINUE 
38 CONTINUE 

CALCULATION OF CFHIC? 

CFFLAG(5*l>-0.0 
DO 2001 N-liS 

IFICFFLAGINI. LE.0.0 1 CFCORRINI-O.O 
IF(CFFLAGIN).lE.O.O) GO TO 2000 
I FICFFLAGIN + l ) .LE.0.0 I GO TO 1000 

CFCORRI Ml BREDEL2I N)*l I ALOGI RE0EL2IMI I-AL0GIREDEL2IH41 1 1 l/( NEXUGGIM 
1 l-REXUGGtN*ll ) ) *1 REXUGGI HI*KHIH 1*1.01 
GO TO 2000 

1000 IFIM.EQ.il CFCORR ( M I »0» 0 
IFIH.EQ.il GO TO 2000 
IFICFFLAGIM-l I. LE.0.0 1 CFCORAIN I “0.0 
IFICFFLAGIM-ll. LE.0.0 I GO TO 2000 

CFCORRI N)*REDEL2(M)*I I AL0GIRECEL2IM-1 ) I-AL0GIRE0EL2I HI 1 1 /(REXUGG f M 
l-D-REXUGGINI l)*<REXUGGIMI*KM(Ht»1.0l 

2000 CONTINUE 

2001 CFNICR(MI-CF2NMMI*CFC0RR(NI 

Cl DETERMINATION OF CORRESPONDING VALUES OF U4, Y + , P», A NO VO* 

DO 3001 M-l.S 

VPLUS1 IMI ■FMIM) /SQRTICF2SLIN 1 1 

PPLUSUN1— 1.0*KM(MI/(CF2SUM)*SORTICF2SUMII) 

VPLUS2INI-0.0 
PPLUS2I r'D'0.0 

IF ICF2N1 1 Ml .Cl .0.0 I VPLUS2I HlnFNIHI /SQRT(CFMECRIN) I 
IFICF2MIINI.GT.0.0I PPLUS2 IN)«-1.0«KNIM)/ ICFMICRIH1 *SQRT(CFHICRI N) 
111 

OP-IIMI 

00 Pi a-l.oP 

UPLuSlI NtR I ■UUGIM .R l/S0RT(CF2SL( Nl I 

YFLUS1 CMiR l«Y I M.R)*REXUGG(N 1*SQRT(CFZSL(NI I All M I 

UPLUS2I N.RI *0.0 

YFLU$2(M,RI*0.0 

IF ICF2HI (M I.GT .0.0 I UPLUS2IN.R l-UUGIM.RI ASORTICFMICRt HI I 
IFICF2MK NI.GT.0.0 1 YPLUS2IH.RI -YIM.R l*REXUGG(N l*SQRT(CFMICR(MI I / 
1XIMI 

37 CONTINUE 
3001 CONTINUE 
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c 

C V11U s » I NT IK ; 'TF riLCUL^TEC BU'VniE' 

c 

C A. SETUP 3AT* 

C 

MRITEI6.7I 

56 FCR^ATUM , 12HRUN NUM8FR! .16, In-, 111 
W«lTtt*i,56l D/TF.RUNNO 

ST FORMAT!//,!] *,6HTGAS» ,F5.2,2F " . AM .AmRHUM- iFA.?, 8M, 19HN0CT* L»M/ 
1ISEC-FT2),//1 
w«ITE(6.57l Tr,,«HU* 

sj fcr‘HT(tth i mil tcavg 1 1 1 util rmo«tn Kill 

I MD3TII! 'Mil 

11*16,531 

I JK-1 

FA FORNATHM,3HlN.,SM,7HtDEG F I , At ,6HF T/SEC , AM , TM 1/SEC ,/ 1 
WRITE (6.A91 

si FORMAT! IX, I*, 3M,F5.2,AM,F5.2 ,6M , F6. 2 , A M ,F6. 3 , 3X.Ein.lt 3X.FT.A, St ,e 
IB. 5) 

HR IT* 16,59 | IJK.MSMJM ,T AVR ( I JK 1 , UGI 1JK I ,DUOK 1 1 JK 1 ,KI 1 JK I * 

IMOOTl tJKI ,Fl I JO 
DO 81 1=1,23 
I J*2*l 
1 J**1 J»1 
LK»I I JK,l I /2 

60 FORMAT! >M,I2,3M.-5.2,1SM,F6.2,A«»F6.J,3M,E1*;.3I 
UR 1 TE I 6 ■ * 3 I IJ.MSI I Jl ,UGt U1 ,CUOMI I JI.KM J1 

WRITE (6,59 I ! JK, MS 1 1 JK I .TAVGt LK I ,U(j t UK I , DUO Ml IJKI.Kt I JK t , HOOT t L K 

II ,FILKI 
fll CONTINUE 

UR I TE 16,31 
C 


r 

B m SUGARY OF PROFILE PARAHETFRS 




t 

51 FORMAT I /////////* 2<tl7H RIH 

NUH6 E R ! 

,i6,iM-,m 



WRlTEf6,51J ^TF ( RUMND 





52 F CRMAT 1 / / 1 1 12H H * 

HREX 

UGG K 

F 


1 CELT* NREDELT? 14 

OEL 

CF2I5LI CF2IMI.CI1 



uRJTEffr«52J 





5 3 FORMAT! 10 M,3Ht I. , 1 TM , 6HF T/SEC ,2 2 M , 3H! N. , 2SK , 3H I N. , / 1 
WRITEI6.53I 
00 55 *=!,* 


5A FORMAT ISM,, , , . F 6. 3 , 3 X , E 1. 3 ,7 K , F6. 2 , 1 M ,£ 10 . 3 , 1 X ,F R. 5 ,2 X , F6 , A , AM , 

IF6.1, 3M.FS. j,4.,F5.3,3X,F7.5,3M,F7.5[ 

UR I TE ! 6, 5A I M,M!M> , REXUGGt H > ,1>GG( N > , KM! Ml ,f HI H ) , OEL T? (M ) , REDEL 2 1 Ml 
1 , H (Ml , DEL ( M I ,C F2SL 1M 1 ,CFM ICR I M I 
55 CONTINUE ' 

UR I TE 16 ,31 

UR ITEI6 ,51 ) DATE, RUN NO 

1005 FORMAT!//, 1UH H M NREX K NREDELT2 

1 CF2I SL ,11 CF21*L,21 C F2 ( MI , A I CORfi. CF2(M1,CI1 

URITEI6, 10051 

1006 FORMAT! 13H IN. ,/l 

WRITE (6, 10061 

l COT FORM ATI AM, [ 2 , 3X ,F 6. 3 ,2 M ,E9. 3 ,2 M , E10 . 3, AM, FT.l , 5M ,F7.5 ,5X,FT. 5, 6M , 
1FT.5,AM,FB.5,5M,FT. 51 
00 100B F = 1 , S 

WRITE 16,10071 K,X ( Ml , REMUGGI Ml ,KMI Ml ,RE0EL2!MI .CF2SLIHI .CF2SL2IM I, 
ICF2M t ( M I tCFCORR (M I , CFM1CR (M I 
1008 CONTINUE 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

E 

[ 

1 

E 

L 

I 

I 


L Ct INDIVIDUAL PROFILE DATA 
C 

DO TC 4*1.1 
*RITE(».3) 

WRI Tt ( 6 .56 I OxTE.RLNNO 
T3 FORMAT! 3H I 
WR1TE16.T3I 

62 FORN4TISX.3HW* ,tl,3*,3HX* ,F6.3,4H I N. t 3X.SH>1GG* .F6.2.TH FTAS'C. 
18X.7HDELT2* .F6.4.4H IN. t 3X, 9tC F 2 1 SL I ■ , F 7, 5 . 3 * , 11 HPPLUS I St I ■ ,F(J. 
25) 

WRITE (6.62 I M,X(N| ,UGGI M) ,OFLT2(H) ,CF2$L(M) .PPLUSIIM) 

63 FCHMAT<12(.3HK> , E1C.3.3X . 6H4COT* ,F7.6tl3K LBN/FTT-SEC .10HWEOEL 
1T2* ,F6.1t6F.RHCF2(WI)> .FT. 5 1 5 * , 1 5HVT EROPLUS ( SL )■ .FB.SI 

WRITE (6, 63 I KWI HJ , MDOT>*(« I ,R:CFL2<N1 ,CF2N1 1*0 , VPLUS1 < Ml 
66 FQRNATI12K.3HF* .FB.5.5X. THVIt RO* .fa.S.TH FTASEC .4X.3HH- ,F5. 3.12 
1X.1CHCF2IMIC I- ,F7.5,6«,13MFP11JS(MI,C1- tfB.M 
WRITE! 6.66 I FMM|,V26KDHIMI,HIMI,CFMICR(M),PPLUS2(M| 

65 FORMAT! 20X.6HNREX* , E9.3 » 11 * .5HQEL- . F3. 3 , 31 X . 1 7HVZER0PLUS I MI.C I ■ 
1.FB.5.//I 

WRITE (6 ,651 REXUGGI * I.DEL (M I, VPLUS2IMI 

66 FORMAT! 5*. 110H RYU YOEL UUG VPLUS 

1ISL) UPLUS($LI VPLUS !M II UPIUSMII SUNl SUM2I 

WRITE 16.66) 

671 FORMAT! 1 3X» SHIN. .AX. 6HFT /SEC . 51 X , THI CORR. I .4X .THICORR. 1 ,7X , 3HI N. .7 
1X.3HIN. I 
WRITE! 6*671 1 
HR ITEI6.75) 

DP-KKI 
ISKIP-5 
00 69 R*l,DP 

68 FOR*AT!5X,I2,*X,F6.6.3X.F6.2,6X,F6.6, 3X.F6.4.6X.FT.2.4X.F5. 2.8X.FT 
1.2,4X,FS.2,7X,FT.S,3X,F7.5I 

WRITE 16. 68 I R.VIM.RI.UIM.R) . YCELIM.RI .UllGIH.RI . YPLUSHM, R I , UPLU* 1 
1 (M.R I ■ VPLU$2(N.RI .UPLU52! M.R! tSUMl IM.R I .SUM2IM.R I 
IFIR.LT.tSKIPI GO TO 69 
WRITE (6. 751 
I $KIP*I jKl P* 5 

69 CONTINUE 

70 CONTINUE 
WRITE! 6.31 

C 

C IX. PUNCHING CF CALCULATED QUANTITIES 
C 

200 FORMAT II6ilX.Il.lX. II. 1X.F6. 3.E10.3 , 1X.F8.5. lX t F6.2.1X,F7.6,1 X.F 

18.3ax.E9.3l 

206 FCRMATI F6.4aX.F6.lt 1X.F5.3. IX. F9. 3) 

205 F0RMATlF7.5ax.FS.5aXtFB.5ax,F7.5aK.F8.Sax.P8.5l 
2011 FORMAT! F7, 5 .1X.F8. 5.1 X.F 7. 5.1 X.F 3.0) 

201 FORMAT I I2.1X.F6. 4. IX. F6.2 . lX.F6.4.1X.F6.4ax,FT. 2.1X.F5.2.1 X.F 7, 2. 
11 X.F 5. 2 .1X.FT.5 i1X.FT.5I 

00 203 M-l.S 

PUNCH 230. DATE .RUNNO.M.X INI , KM! Mt . FMIMl .UGGIN I .NOOTMIM t , VtEROMI Ml 
1.REXUGGIM1 

PUNCH 206, DELTZtM I .REDEL2IN I .HIM) .DEL I Ml 

PUNCH 203. CF2SL!NI .PPLUS1INI .YPLUS1 INI .CF2MI (Ml .PPLUS2IM) . VPLUS2! 
INI 

1FICF2MI !M I.EQ.0.0 ISFlAGI H l*0.0 
IFICF2M1 (MI.NE.O.OI SFLAGIMI-1.0 

PUNCH 2011. CF2SL21NI.CFC0RRIN) ,CFMICR!MI .SFLAG(N) 

0P-2INI 
00 202 R-l.OP 

PUNCH 201. R.VIM.RI.UIM.R I. YOEL (M.R I .UUG! M.R) , YPLU511M.R I .UPLUS1 (M 
1, Rl. YPLUS2 I M.R) , UPLU52 (N, R I, SUM1I M.R I, SUN2 IN. R I 

202 CONTINUE 

203 CONTINUE 
C 

lFIRUNN0.LT.3t GO TO 2 

RETURN 

END 

A DAT A 
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APPENDIX C 

SUMMARY OF SEMI-EMPIRICAL REPRESENTATIONS OF DATA 


Assumed mixing length distributions: 

i. Two- layer model 



y* < y+ 

O 

II 


y c /5 < y/B < a/k 

>> 

li 

(k = 0.44) 

y/5 > X/k 

i = X5 


li. Continuous modified Van Driest 

model 

y/6 < X/k 


/y + yr* r \ 

>> 

IJ 

1 ■ exp \~^r) 


1 

L— _J 

y/d > X/k 

to 

II 


Truncation of mixing 

length: 



- 1/8 

X = 0.25 Re &2 (1-67,5 F) 

where X is truncated at a value of 0.085* 

c) Shear stress distribution assumed in inner regions of 
layer : 


t + = 1 + uV + p + y + 

1 . i 
y 

f(* 

X 2 

:) dy 


- 

J o 



] 
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d) Correlation of = y* J t + and A. with and p + : 


^ * J c| p+=0 U - V^’ 


A, - A, I [1 - <U<)P + ] 

|p+=0 

where 



11.0 - 18.0 V + 
w 


and 


|p + =o 


the functions 


26.0 - 88.0 v+ + 110.0 (V+) 2 

W rl 

26.0 - 88.0 v; + 210.0 (v*) 2 

Qy (V^) and Q*(V*) are given 



> 0 



< 0 


by 


- 0.08 

0.0 

0.005 

0.01 

0.015 

0.02 

0.025 

0.03 

0.035 

0.04 

0.045 

0.05 

0.06 

0.4 


M V w> 


63.3 

19.7 

18.0 
16.6 

15.7 

15.0 

14.5 

14.1 

13.8 

13.6 

13.4 
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